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A new experimental investigation of decaying turbulence generated by a low-blockage 
space-filling fractal square grid is presented. We find agreement with previous works by 
Seoud & Vassilicos [ "Dissipation and decay of fractal-generated turbulence" , Phys. Flu- 
ids 19, 035103 (2007)] and Mazellier & Vassilicos ["Turbulence without the Richardson- 
Kolmogorov cascade", Phys. Fluids 22, 075101 (2010)] but also extend the length of 
the assessed decay region and consolidate the results by repeating the experiments with 
different probes of increased spatial resolution. It is confirmed that this moderately high 
Reynolds number Rex turbulence (up to Re\ ^ 350 here) does not follow the classi- 
cal high Reynolds number scaling of the dissipation rate e ~ u' 3 /L and does not obey 
the equivalent proportionality between the Taylor-based Reynolds number Re\ and the 
ratio of integral scale L to Taylor micro-scale A. Instead we observe an approximate pro- 
portionality between L and A during decay. This non-classical behaviour is investigated 
by studying how the energy spectra evolve during decay and examining how well they 
can be described by self-preserving single-length scale forms. A detailed study of homo- 
geneity and isotropy is also presented which reveals the presence of transverse energy 
transport and pressure transport in the part of the turbulence decay region where we 
take data (even though previous studies found mean flow and turbulence intensity pro- 
files to be approximately homogeneous in much of the decay region). The exceptionally 
fast turbulence decay observed in the part of the decay region where we take data is 
consistent with the non-classical behaviour of the dissipation rate. Measurements with a 
regular square mesh grid as well as comparisons with active grid experiments by Myd- 
larski & Warhaft ["On the onset of high- Reynolds- number grid-generated wind tunnel 
turbulence", J. Fluid Mech. vol. 320 (1996)] and Kang, Chester & Meveneau ["Decaying 
turbulence in an active-grid-generated flow and comparisons with large-eddy simulation" , 
J. Fluid Mech. vol. 480 (2003)] are also presented to highlight the similarities and differ- 
ences between these turbulent flows and the turbulence generated by our fractal square 
grid. 



1. Introduction 

At high enough Reynolds numbers, the local viscous dissipation rate e of the lo- 
cal average turbulent kinetic energy K scales with K and a local correlation length 
scale L, i.e. e ~ K 3 / 2 /L . At least, this is what one reads in turbulence textbooks (see, 



for example , iBatchelor 1953L Tennekes k, Lumley r 1972 Lum evl 1992 \_ Townsendl 1956 



Frischlll995t lLesieur)ll997t|Mathieu fc Scotdl200d: |Popell200d: ISagaut fc Cambonll2008l ) 
Tennekes fc Lumlevl (|1972l ) introduce this scaling in their very first chapter with the 



Townsendl dl956ft 



words "it is one of the cornerstone assumptions of turbulence theory" 
uses it explicitly in his treatment of free turbulent shear flows (see page 197 in I Townsendl 
19561) which includes wakes, jets, shear layers, etc. Since G.I Taylor introduced it in 
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1935 (|Tavlorlll935| ). this scaling i s also customar i ly used in theories of de caying homo- 
geneous isotropic turbulence (see Batchelorlll953 : Frischl[l995l Rotta 1972) and in anal- 
yses of wind tunnel simulatio ns of such turbulence (e.g. Batchelor fc Townsendl 1948 
Comte-Bellot fc Corrsinlll966l ) in the form 



L 



(1.1) 



where u' is the r.m.s. velocity fluctuation, L is an integral length scale and C £ is a con- 
stant independent of time, space and R eynolds number when the Reynolds number is 
large enough. However, as iTavlorl (Il935l) was careful to note, the constant C e does not 
need to be the same irr espective of the boundaries (initial conditions) where the tur- 
bulence is produced (se e Burattini. Lavoie fc Antonia 2005 : Mazellier fc Vassilicos! 2008 : 
Goto fc Vassilicosll2009h . 

In high Reynolds number self-preserving free turbulent shear flows, the cornerstone 
scaling e ~ K 3 / 2 /L determines the entire dependence of e on the streamw ise coordinate 
and ascertains its independence on Reynolds number (see lTownsendl[l95l . This corner- 
stone scaling is also effectively used in turbulence models suc h as K-e (see lPopdEoobT ) 
and in Large Eddy Simulations (see Lesieurlll997t |Popell2000l ) . The assumption that e is 
independent of Reynolds number when the Reynolds number is large enough is an in- 
separ able part of the Richardson-Kolmogorov cascade ( Tennekes fc Lumlev 19721 : iFrischl 
1995). This is the celebrated nonlinear dissipation mechanism of the turbulence whereby, 
within a finite time Lj\/~K (the same time scale for all high enough Reynolds num- 
bers), smaller and smaller "eddies" are generated till eddies so small are formed which 
can very quickly lose their kinetic energy by linear viscous dissipation. The higher the 
Reynolds number, the smaller the size of these necessary dissipative eddies but the time 
scale L/\/K for energy to cascade to them from the large eddies remains the same. The 
dissipation rate e is proportional to K divided by this time, and therefore e ~ K 3 / 2 /L. 

In various high Reynolds number self-preserving free turbulent shear flows as in wind 
tunnel grid-generated turbulence, K and L vary with streamwise downstream distance 
x-xo (where xo is an effective/ virtual origin) as power laws. Specifically, K ~ U^ ( x L ^" )~ n 
and L ~ Lsi^j^ 3 -)" 1 where Lb is a length-scale characterising the inlet and Uoo is the ap- 
propriate inlet velocity scale. In table [1] we recall the generally accepted values taken by 
the exponents n and in in plane wakes, axisymmctric wakes, self-propelled plane wakes, 
self-propelled axisymmctric wakes, mi xing layers, plane jets, axisymmetric jets and wind- 
tunn el grid- generated turbulence ffrom lTennekes fc Lumlevlll972tlComte-Bellot fc Corrsin 
19661) . Estimating a Taylor microscale A from e ~ vK/X 2 where v is the kinematic vis- 
cosity of the fluid, and then applying the cornerstone assumption e ~ K 3 ! 2 /L to all these 
flows yields the following two relations: 



■''o 



(1.2) 



and 



xo )f . 



(1.3) 



where Reo = UooLb/v is the inlet Reynolds number and yKX/v is a local Taylor 
microscale-based Reynolds number. The different values of y — ^ are given in table [TJ 
Remarkably, e ~ K 3 / 2 /L implies that L/X~ \f~KX\v in all these flows whatever the values 
of n and m, meaning that L/X ~ \f~KX\v collapses the Reo and the x dependencies in the 
same way for all these flows. We stress that this collapse is the immediate consequence 



Decay of multi-scale grid turbulence 



3 



Table 1. Powers law exponents characteristic of the downstream evolutions of K, L, L/X 



K L L/X 



Plane wake 1 1/2 

Axisymmetric wake 4/3 1/3 -1/6 

Self-propelled plane wake 3/2 1/4 -1/4 

Self-propelled axisymmetric wake 8/5 1/5 -3/10 

Mixing layer 1 1/2 

Plane jet 1 1 1/4 

Axisymmetric jet 2 1 

Regular grid turbulence 1.25 0.35 -0.14 



of e ~ K z l 2 /L. The relation L/X ~ \f~KX/v simply reflects the Richardson-Kolmogorov 
cascade: the higher the Reynolds number, the smaller the size of the dissipative eddies, 
i.e. the greater the range of e xcited scales and the greater L/X. 

As noted bv lLumlevI dl992h . by 1992 there had not been too much detailed and com- 
prehensive questioning of data to establish the validity of e ~ K 3 / 2 /L but he wrote: "I 
hardly think the m atter is really much in question". He cited the data compilations of 
Sreenivasanl (| 1984T ) which suggested that C e does become constant at Re\ = u'X/v larger 
than about 50 for wind tunnel turbulence generated by various biplane square-mesh 
grids, but there seemed to be little else at the time. Since then, direct numerical sim- 
ulations (DNS) of high Reynolds number statistically stationary homogeneous isotropic 
turbulence have significantly strengthened s upport for the constanc y of C E at Rex greater 
than a bout 150 (see compilation of data in lBurattini et al\ ( 20051 ). see also Sreenivasanl 
( 19981 ) ) . Other turbulent flows have also been tried in the past fifteen years or so such as 
vari ous turbulent wakes and jets and wind tunnel turb ulence generated by active grids 
(see Burattini et al. 2005 : Mazellier fc Vassilicosl [20081) with some, perhaps less clear, 
support of the constancy of C F at large enough Re \ (perhaps larger than about 200 if L 
is defined appropriately, see iBurattini et al\ (|2005l )) and also some clear indications that 
the high Rey nolds number constant value of C e is not universal, as indeed cautioned by 
Tavlorl (|l935h . 

A decade ago, iQueiros-Conde fe Vassilicos ( 200ll ) took the opposite approach and 
asked whether it might be possible to break e ~ K 3 I 2 /L in some fundamental way 
in some flows, and so they proposed generating turbulence with fractal/ multiscale ob- 
jects/stirrers/inlet conditions. Some years later. iHurst fc Vassilicosl (|2007t published an 
exploratory study of wind tunnel grid-generated turbulence where they tried twenty one 
different planar grids from three different families of passive fractal/multiscale grids: frac- 
tal cross grids, fractal I grids and fractal square grids. They ascertained that the fractal 
dimension Df of these grids needs to take the maximal value Dt = 2 for least downstream 
turbulence inhomogeneity. They also identified some important grid-defining parameters 
(such as the thickness ratio t r , see figure [1] and table [3|) and some of their effects on the 
flow, in particular on the Reynolds number Re\ which they showed can reach high values 
with some of these grids in small and conventional sized wind tunnels, comparable to 
values of Rex achieved with active grids in similar wind tunnels and wind speeds. Their 
most interesting, and in fact intriguing, results were for their space- filling (Df = 2) low- 
blockage (25%) fractal square grids (see figure[T]). Fractal square grids have therefore been 
the multisca le grids of choice in most subsequent works on multiscale/fractal- generated 



turbulence (Seoud fc Vassilicos 20071 Nagata. Suzuki. Sakai. Havase fc Kub 



l-gcncratcd 
olboOSftO: 



Stresing. Peinke. Seoud fc Vassilicosl2010l : lMazellier fc Vassilicosl2010USuzuki. Nagata. Sakai fc Ukai 
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2010t Laizet fc Vassilicos] 2011 ). For the case of space- filling low-blockage fractal square 
grids , iHurst fc Vassilicosl (|2007t ) found a protracted region between the grid and a dis- 
tance Xp e ak downstream of the grid where the turbulence progressively builds up; and 
a decay region region at x > x pea k where the turbulence continuously decays down- 
stream. They reported a very fast turbulence decay which they fitted with an expo- 
nential and also reported very slow downstream growths of the l ongitudinal and lateral 
integr al length-scales and of the Taylor microscale. (Very recently, iKrogstad fe Davidson! 
( 20111 ) studied the decay behind multiscale cross grids and found conventional decay 
rates. Note that for multisc ale cross grids our prior publications did not claim fast, 
unconventional, decay rates (jHurst fc Vassilicosl 120071) . This may serve as further jus- 
tification for focusing attention on fractal square grids in the present paper. Even so, 
multiscale cross grids ha ye been used successfully in some recent studies for enhan cing 
the Reynolds number, seelKinze . Wolf. Holzner. Luthi. Tropea fc Kinzelbach ( 2011 ) and 
Geipel Henry Goh fc Lindstedtl (|2010l) .) 

Seoud fc Vassilicosl ( 2007 ) concentrated their attention on the decay region of turbu- 



le nce generated by space- fill ing low-blockage fractal square grids and confirmed the results 
of Hurst fc Vassilicosl ( 20071 ). In particular, they showed that L/X remains approximately 
constant whilst Re\ decays with downstream distance x and they noted that this be- 
haviour implies a fundamental break from (jl.ip where C £ is constant. They also found 
that one-dimensional longitudinal energy spectra at different downstream centreline lo- 
cations x can be made to collapse with u' and a single length-scale, as opposed to the 
two length-scales (L and Kolmogorov microscale) required by Richardson-Kolmogorov 
phenomenology. Finally, they also carried out homogeneity assessments in terms of vari- 
ous profiles (mean flow, turbulence intensity, turbulence production rate) as well as some 
isotropy assessments. 

Mazellier fc Vassilicos ( 2010l ) also worked on wind tunnel turbulence generated by 



space-filling low-blockage fractal square grids. They introduced the wake-interaction 
length-scale x* which is defined in terms of the largest length and thickness on the grid 
and they showed from their data that x pea k « 0.5x*. They documented how very inhomo- 
geneous and non-Gaussian the turbulent velocity statistics are in the production region 
near the grid and how homogeneous and Gaussian they appea r by comparison beyond 
0.5a;*. They confirmed the findings of Hurst fc Vassilicosl ( 20071 ) and lSeoud fc Vassilicosl 



( 20071 ) and added the observation that both Re\ and L/X are increasing functions of the 
inlet velocity Uoo- Thus, the value of L/X seems to be set by the inlet Reynolds number, 
in this case defined as Reo = UooX^/v for example. 

Finally, iMazellier fc Vassilicosl (|2010| ) brought the two differe nt single-scale turbulence 
decay behaviours of George! (Il992l ) and lGeoree fc Wanl (l2009h into a single framework 
which they used to analyse the turbulence decay in the downstream region beyond x pea k « 
0.5x*. This allowed them to introduce and confirm against t heir data the notions that , 
in t he decay region, th e fast turbulence decay observed bv IHurst fc Vassilicosl ( 2007 ) 
and ISeoud fc Vassilicosl ([20071 ) may not be exponential but a fast decaying power-law 
and that L and A are in fact increasing functions of x which keep L/X approximately 
constant. 

The results of lHurst fc Vassilicosl (|2007l) . ISeoud fc Vassilicosl (|2007l) and lMazellier 



( 2010l ) suggest that, in the decay region downstream of space-filling low- blockage fractal 
square grids, high Reynolds number turbulence is such that 



Vassilicos 



L/\~Re%A(?—^) 



(1.4) 
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and 



Re x ~Re%B('- 



x . 



(1.5) 



where A is a slow- varying dimcnsionless function of ^r 2 - (in fact effectively constant), B 
is a fast-decreasing dimcnsionless function of (perhaps even as fast as exponential) , 
and a and /3 are positive real numbers. 

Assuming that the dissipation-scale tu rbulenc e struc ture is approximately isotropic, we 
now use the relation e = I5vu' 2 /X 2 which iTavlorl (|1935l ) obtained for isotropic turbulence. 
With ([l.ip this relation implies 

and, clearly, C £ cannot be constant (independent of Reo and x) with Reo and x depen- 
dencies of L/X and Re\ such as those observed in wind tunnel turbulence generated by 
space-filling low-blockage fractal square grids. Instead, 



C E = WRe 



~^A( X - — 



(1.7) 



which means that C e should be increasing fast in the downstream direction but which 
also means that a plot of C £ versus Re\ can be quite different depending on whether 
Re\ is varied by varying Reo whilst staying at the same position x or by moving along x 
whilst keeping Reg constant. This is a point which we discuss and attempt to bring out 
clearly in the present paper. 

Relations (|1 .4[) and (|1.5|) and their conseq uent decoupling of L/X an d R e\ were ob- 

served at moderate to high values of Re\ where lSeoud fe Vassilicosl (|2007h and lMazellier fe Vassilico'sl 

also observed a well-defined broad power-law energy spectrum. Indeed Re\ needs 
to be lar ge enough for t he study of fully developed turbulence. Active grids were intro- 
duced by iMakital ( 199lh to improve on the Reynolds number values achieved by regular 
grids in conventional wind tunnels. Fractal square grids achieve comparably high values 
of Re\ but also a far wider range of Re\ values along the streamwise direction. This 
makes if much easier to study i?e,\,-depcndcncics, a point which we make and discuss in 
some detail in the present paper. 

In this paper we report an experimental assessment of turbulent flows generated by 
a low-blockage space-filling fractal square grid (see figure [1]) and a regular square-mesh 
grid. The main focus of this paper is to complement former research on fractal-generated 
turbulence by extending the assessed decay region and using the new data to re-address 
the previously reported dramatic departur e from C e = Const and A(^ 2 -) = B(^r 2 -) 
and the abnormally high decay ex ponents ( Hurst fe Vassilicos 2007t ISeoud fe Vassilicosl 
2007t IMazellier fe Vassilicosl [ioToh . We provide estimates of these exponents, and also 
show that a « (5 and that our fractal-generated turbulen ce behaves in a way which i s 
very close to self-preserving single-length scale turbulence (|Mazellier fc Vassilicosll2010l) . 



particularly if the turbulence anisotropy is taken into account when calculating 3D energy 
spe ctra. We also show that , even though previous studies bv lSeoud fc Vassilicosl (|2007t ) 



and lMazellier fc Vassilicosl ( 201dl ) found that the mean flow and turbulence profiles are 
approximately homogeneous in much of the decay region, there nevertheless remains sig- 
nificant transverse turbulent transport of turbulent kinetic energy and turbulent trans- 
port of pressure. The decaying turbulence is therefore not homogeneous and isotropic in 
terms of third order one point statistics even though it more closely is in terms of lower 
order one point statistics. Whenever possible a comparison between fractal-generated 
and non-fractal-generated turbulence is made emphasising similarities and differences. 
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In the following section we describe the experimental apparatus as well as the anemom- 
etry systems, probes and the details of the data acquisition. The experimental results 
are presented in Sec. [3] and are organised in four subse ctions. In Sec. 13.11 it is sugges ted 
that the wake-interaction length-scale introduced by iMazellier fc Vassilicosl (|2010h to 
characterise the extent of the production region in the lee of the fractal grid is also 
meaningful for regular static grids. In Sect. 13.21 the homogeneity and isotropy of the 
frac tal-generated flow is inve s tigate d following the methodology used by ICorrsinl (|1963l) 
and IComte-Bellot fe Corrsinl (| 1966T) for regular static grids. In Sec. 13.31 & 13.41 the nor- 
malised energy dissipation rate and the decay law are re-assessed using the new data. In 
Sec. 13.51 we investigate the possibility of a self-similar, single-length-scale behaviour by 
collapsing the ID energy spectra and the 2™ d -order structure functions using large-scale 
variables; also the 3D energy spectrum function is calculated to provide isotropy correc- 
tions on the collapse. In Sec. |4]we end this paper by highlighting the main conclusions 
drawn from the present measurements and discuss some of the questions raised. 



2. The experimental setup 

2.1. Experimental hardware 

T he experiments are per f ormed in the T = 0.46m wind-tunnel described in some detail 
in Mazellier fe Vassilicosl (|2010h and sketched in figure [5] (T is the lateral width of the 



tunnel's square test section). The inlet velocity Uoo is imposed and stabilised with a PID 
feedback controller using the static pressure difference across the 8:1 contraction and 
the temperature near the inlet of the test section which are measured using a Furness 
Controls micromanometer FCO510. 

All data are taken with one- and two-component hot-wire anemometers operating in 
constant-temperature mode (CTA). The hot-wires are driven by a DANTEC StrcamLinc 
CTA system with an in-built signal conditioner. We use both square- and sine-wave 
testing to measure the cut-off frequency at the verge of attenuation (f cut-off) anc ^ at the 
standard '-3dB' attenuation level (f cutoff)- ^ n table [2] we present the results from the 
electronic testing of our anemometry system. Further information concerning electronic 
testing of thermal anemometers an d a discussion of the consistency between the square 
and sine- wave tests can be found in iFrevmuthl (|1977l ) . 

For the single component measurements three different single- wires (SW) are used with 
a sensing length (l w ) of lmm, 0.45mm & 0.2mm respectively. For the two component 
measurements two cross- wires (XW) with sensing lengths of l w = 0.5 & lmm respectively 
are used, but for both the separation between the wires is around lmm. All the sensors 
except the l w = lmm XW are based on Dantec probes modified to use in-housc etched 
Platinum-(10%)Rhodium Wollaston wires soldered to the prongs (further details can be 
found in tabled]). The l w = lmm XW is a Dantec 55P51 tungsten probe. It should be 
noted that the l w = 0.2mm single-wire, which has a diameter of d w = 1.27/im, is oper- 
ated in the limit of the bridge stability, on the verge of having non-damped oscillations. 
Nonetheless, the sine-wave test indicated that /°u t B off was about 40kHz. The hot-wires 
are calibrated at the beginning and at the end of each measurement campaign using a 
4* ft, -order polynomial in the SW case and a velocity-pitch map in the XW case. Note that, 
unless otherwise stated, the data shown are acquired with the l w « lmm SW hot-wire 
probe. All the two-component data presented are acquired with the l w « 0.5mm XW 
except the spanwise traverse data presented in Sec. 13.2.11 

Note that two other anemometry systems have been used as well in order to compare 
with previous experimental results, but these results are not included here. The other 
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SW/XW l w d w l w /d w Hot-wire U„ /^ fl / c ° u d t s off l w /r, 
(mm) i^m) probe (rns -1 ) {kHz) (kHz) 



SW ~ 1 5.1 196 55P16 

SW - 0.45 2.5 180 55P16 
SW - 0.2 1.27 157 55P11 



10 - 25 - 12 7-3 

15 - 32 - 16 9-5 

10 -45 -21 3-2 

15 - 45 - 23 4-2 

10 > 50 - 40 - 1 

15 > 50 - 40 2-1 



XW - 0.5 2.5 200 55P51 
XW 1.0 5 200 55P51 15 -30 



10 -45 -21 3-2 

15 - 45 - 23 4-2 



Table 2. Details on the hot-wires, cut-off frequencies & resolution. l w and d w are the sensing 
length and diameter of the wires, l w /n is the ratio between the sensing length and the Kolmogorov 
inner length-scale and Uoa is the inlet velocity. /cut_ fi i s the cut-off frequency corresponding to 
-3dB signal attenuation and /cut-off i s the highest frequency with negligible attenuation. 



anemo metr y systems are: the AALab AN-1005 CTA system used in Seoud fc Vassilicosl 
( 20071) and kurst fc Vassilicosl (l2007t) and the PISA 55M10 C TA bridge with a DISA 
55D26 signal conditioner used in Mazellier &: Vassilicos ( 201dt ). It is found that the re- 
sults obtained with the DISA 55M10 CTA unit closely match those obtained with the 
StreamLine CTA system, when the same hot-wire probe is used, except at very high 
frequencies where the higher noise floor of the DISA CTA system buries the velocity 
signal. On the other hand it is found that the measurements taken with the AALab 
AN-1005 CTA system are significantly different at frequencies above 6kHz and therefore 
the turbulence statistics involving velocity derivatives are significantly different. This is 
likely the reason for the difference be t ween the normalised energy dissipation rate C e 
results reported in lSeoud fc Vassilicosl (|2007l ) and the ones presented in Sec. 13.31 of this 



paper. The comparison between the results of the different anemometry systems will be 
presented elsewhere. 



2.2. Data acquisition and signal processing 

The pressure and temperature measurements are digitally transferred to the computer 
using a parallel port. The analogue signal from the anemometers is sampled using a 16- 
Bit National Instruments NI-6229(USB) card, at a sampling frequency set to be higher 
than twice the analogue low-pass filtering frequency (30kHz). The data acquisition and 
signal processing are performed with the commercial software MATLAB™ . 

The turbulent velocity signal was acquired for 9min corresponding to more than 
100,000 integral-time scales. This was confirmed to be sufficient for converged measured 
statistics of interest such as the integral scale, the first four moments of the velocity 
signal and the 2 nd moment of the velocity derivative signal. The time-varying turbu- 
lent signal was converted into spati ally- varying by m e ans o f a local Taylor's hypothesis 
following the algorithm proposed in iKahalerras et all (]1998f ). Before Taylor's hypothesis 
is used the signal is digitally filtered at a frequency corresponding to k\r\ ~ 1.1 (where 
■q = (i^/e) 1 / 4 is the Kolmogorov inner length-scale and k\ the wavenumber) using a 
4* ft, -order Butterworth filter to eliminate higher frequency noise. 

The integral scale L u is estimated as 



L h 



/ f(r)dr, 
Jo 
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Grid N L to L r t r R L Rt o M eff 

(mm) (mm) (mm) 

SFG 4 237.8 19.2 8 17 0.5 2.57 0.25 26.2 

RG 1 60 10 1 1 1 1 0.32 60 

Table 3. Details of the regular grid (RG) and the low-blockage space-filling fractal square grid 
(SFG). N is the number of fractal iterations of the grids (for a regular grid N = 1). Lq is the 
centreline distance separating the largest bars and io their lateral thickness, see figure[T] L r and 
t r are respectively the length and thickness ratio between the largest and smallest bars. Rl and 
Rt are the length and thickness ratio between two consecutive fractal iterations. Rt is related to 
t r and iV via t r = R£ . The blockage ratio a is defin ed as the fraction of area o ccupied by the 
grid and M e // is the effective mesh size as defined in iHurstll ; Vassilicos (2007) which reverts 
to the definition of mesh size for a regular grid. 



where f(r) = u(x)u(x + r)/u(x) 2 is the auto-correlation function of the streamwise ve- 
locity fluctuations for streamwise separations r and rt is maximum integration range 
taken to be about 10 times the integral length scale. It was checked that (i) changing the 
integration limit tl by a factor between 2/3 and 2 has little effect on the numerical value 
of the integral scale and (ii) the choice of r^, if large enough, does not influence the way 
that L u varies with downstream distance. The transverse integral scale is estimated in a 
similar way. The longitudinal and transverse spectra are calculated using an FFT based 
periodogram algorithm using a Hanning window with 50% overlap and window length 
equivalent to at least 180 integral length scales. The dissipation e is estimated from the 
longitudinal wavenumber spectra F\\ as 



rk m ax 



where k m i n and k max are determined by the window length and the sampling frequency 
respectively. To reduce the unavoidable cont amination of no ise at high frequencies (which 



can bias the dissipation estimate) we follow I Antonial (|2003l) and fit an exponential curve 



to the high frequency end of the spectra which we then integrate. We checked that calcu- 
lating the dissipation with and without Antonia's (2003) method changes the dissipation 
by less than 4% in the worst case. 

It might be worth mentioning that the measurements of the fractal grid-generated tur- 
bulence posed a lesser challenge to hot-wire ancmomctry than the regular grid-generated 
turbulence quite simply because the turbulent signal to anemometry noise ratio is higher 
in the former case, but nonetheless the Kolmogorov microscales (which influence the max- 
imum frequency to be measured) for the highest Re\ measurement location (Rex ~ 350 
and Re\ « 150 respectively) are roughly the same (77 « 0.11mm and r\ « 0.13mm respec- 
tively). 

2.3. Turbulence generating grids 

The bulk part of the measurements are performed on turbulence generated by a low- 
blockage space-filling fractal square grid (SFG) with 4 'fractal iterations' and a thickness 
ratio of t r = 17, see figure CD It is one of the grids used in the experimental setup of 
Mazellier fc Vassilicosl (|2010l) where further details of the fractal grids and their design 
can be found. Measurements of turbulence generated by a regular bi-plane grid (RG) 
with a square mesh and composed of square rods are also performed. The summary of 
the relevant grid design parameters is given in table |3] 
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Figure 1. Low-blockage space-filling fractal square grid (SFG). The grid is space-filling because 
the fractal dimension of its delimiting line takes the maximum value of 2 over the range of scales 
on the grid. In the limit of infinite number of fractal iterations the blockage ratio will tend to 
unity, without taking bar thickness into account. However with only four iterations and the bar 
thickness in the figure the grid's blockage ratio is only 25%. 



3. Results 

The turbulent field in the le e of the space-filling fract a l square grids can be conside red 



to have two distinct regions ( Hurst fc Vassilicoa 2007 : Mazellier fe Vassilicos 2010h : a 



production region where the turbulent kinetic energy (on the centreline) is increasing 
and the flow is being homogenised, and a decay region where the energy of the turbulent 
fluctuations are rapidly decreasing and the flow is roughly homogeneous with an isotropy 
factor around u'/v' ~ 1.1 - 1.25, where v! and v' are the longitudinal and transverse r.m.s. 
velocities respectively 

3.1. The wake-interaction length-scale 

Mazellier fc Vassilicos! (|2010l) introduced the wake-interaction length-scale a;* = L^/to (see 
definitions of Lq & to in figure [T] and in the caption of table [3]) to characterise the extent 
of the turbulence production region in the lee of the space-filling fractal square grids. This 
length-scale is based on the largest square of the grid since the wakes it generates are the 
last to interact, although there is a characteristic wake- interaction length-scale for each 
grid iteration (for a schematic of the wake interactions occ urring at different strcamwise 
locations refer to figure 4a in Mazellier fc Vassilicos 2010l ). They then related the wake 
interaction length-scale with the location of the maximum of the turbulence intensity 
along the centreline x pea k , which marks the end of the production region and the start 
of the decay region and found that x pea k/x* « 0.45. Note that this is not the only peak 
in turbulence intensity in the domain nor is it the overall maximum, but it is the last 
peak occurring; furthest downstrea m before decay. This can be seen for example in figure 
9 in iMazellier fc Vassilicosl <|2010T ). where the streamwise variations of the turbulence 
intensity both along the centreline and along an off-centre parallel line are shown. The 
turbulence intensity along this particular off-centre line peaks much closer to the grid 
and at a higher intensity value than the turbulence recorded on the centreline. 

The wake-interaction length-scale can also be defined for a regular gr id, where the mesh 
size an d the bar thickness are now the relevant parameters, = M 2 /tn. |javesh fc Warhaftl 
( 1992h measured the turbulence intensity very near the grid, 1 < x/M < 30 and observed 
two different regions, a highly inhomogeneous region up to x/M rs 3 which is a produc- 
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tion region where the turbulence intensity increases along a streamwise line crossing half 
distance between grid bars and a decay region beyond that. Not e that x/M « 3 corre- 
spond s to x pea k/x* rs 0.55 close to x pea k/x* ~ 0.45 encountered by iMazellier fc Vassilicos! 
( 2010l ) for the fractal square grids. A qualitatively similar conclusion can be drawn from 



the direct numerical simulation of turbulence genera ted by a regular grid presented in 
Ertunc. Ozvilmaz. Lienhart. Durst fc Beronovl ( 2010l ). In their figure 16 one can find the 



development of the turbulent kinetic energy very close to the grid 0.5 < x/M < 10 along 
three straight streamwise lines located, respectively, behind a grid bar, half-distance be- 
tween bars and in-between the other two traverses. It can be seen that the turbulence 
intensity peaks first directly behind the grid bar at x pea k/M » 1 and lastly behind the 
half-distance between grid bars (somewhat equivalent to the centreline in the square frac- 
tal grid) at x pea k/M » 2.5. This latter streamwise location corresponds to x pea k/x* « 0.68, 
once more not far from x pea k/x* w 0.45. Note nonetheless that this simulation was per- 
formed at very low Reynolds numbers, Re\ < 17, so care must be taken in quantitative 
comparisons. 

Note that x pea k/x* appears to be slightly higher for the regular static grids than for 
the fractal square grids. This is likely due not only to the typically low Reynolds numbers 
generated by the regular grids but also to the characteristic production mechanism of 
the fractal square grids, i.e. before the larger wakes interact all the smaller wakes have 
already interacted and generated turbulence that increases the growth rate of the larger 
wakes, thus making them meet closer to the grid and therefore causing a smaller value 

of X p eak fx*. 

The fact that the fractal grid has multiple wake-interaction length-scales, for the 
present fractal square grid ranging from a few centimetres to more than a meter, is pre- 
cisely part of what makes the fractal grid generate turbulence that is qualitatively differ- 
ent from regular grid-generated turbulence. Consequently one could expect that a fractal 
grid designed so that it produces a narrow range or a single dominant wake-interaction 
length-scale, will lead to turbulence that is similar to regular grid-generated turbulence. 
Irlurst fc Vassilicosl ( 2007 ) included in their study the assessment of fractal cross grids, 



which resemble regular grids but with bars of varying thicknesses. The ratios between 
the thickest and the thinnest bars of their fractal cross grids ranged from 2.0 to 3.3, thus 
yielding a narrow span of wake-interaction length-scales. Furthermo re, the wake inter- 
action pattern of the fractal cross grids, as designed and studied in Hurst fc Vassilicosl 



(|2007h . is considerably different from the wake interaction pattern of their fractal square 



grids. In the fractal square grids case, the main interaction events occur when similar 
sized wakes meet, whereas in the fractal cross grids the main interaction events occur 
between adjacent wakes, which may or may not be of similar size. Therefore one could 
expect the results obtained with fractal cross grids, for example the power-law turbu- 
lence decay exponent, not to be very different from the t ypical results found for reg ular 
grid-generated turbulence. In fact, examining figure 10 in Hurst fc Vassilicosl ( 2007 ) one 



can see that the turbulence decays as {x-Xo)~ n with 1 < n < 1.5 for xq « 0, although they 
encounter a general difficulty of finding the appropriate virtual origin. We will return to 
the problem of finding the appropriate virtual origin and the power-law decay exponent 
in Sec. 13.41 where we present different power-law decay fitting methods applied to our 
data. 

3.2. Homogeneity, isotropy and wall interference 

3.2.1. Homogeneity 

Previous experim ental investigations on the tu rbulence generated by space-filling frac- 



tal square grids, e.g. Mazellier fc Vassilicos ( 2010l ). reported that the flow field close to the 
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grid is highly inhomogeneous. It was also observed that during the process of turbulent 
kinetic energy build up the turbulent flow is simultaneously homogenised by turbulent 
diffusion, and by the time it reaches a peak in turbulence intensity (what they considered 
to be the threshold betwe en the production and de cay regions) the flow has smoothed 
out most inhomogeneities. Seoud fc Vassilicosl ( 2007 ) measured the turbulent kinetic en- 



ergy production in various planes perpendicular to the mean flow along the centreline 
and observed that the turbulent production decreases rapidly just after the peak, i.e. 
where 0.45(« x pea k/x*) < x/x± < 0.75 and that the turbulent energy production typically 
re presents less than 30% of the dissipation and never exceeds 20% beyond this region. 
Mazellier fc Vassilicosl (|2010l) compared the characteristic time scales of the mean ve- 



locity gradients (^r) 1 and (f^-) 1 (where U is the streamwise mean velocity and y is a 
coordinate along the horizontal normal to the streamwise direction) with the time scale 
associated with the energy-containing eddies and reached the conclusion that beyond the 
peak the mean gradient time scale is typically one to two orders of magnitude larger. 
Consequently the small-scale turbulence dynamics are not affected by large-scale mean 
flow inhomogeneities. 

H ere we complement the p r eviou s analyses by following the approach of Corrsinl ( 19631) 



and lComte-Bellot fc Corrsinl (|1966l ) and using some of their homogeneity criteria, as they 



did for regular grids. The commonly accepted 'rule-of-thumb' for the regular grids is that 
the turbulent flow can be considered statistically homogeneous in transverse planes for 
x/M > 30 and the unavoidable inho mogeneity alo ng the mean flow direction becomes 
relatively unimportant for x/M > 40 ( Corrsinlll963l ). 



For the downstream decaying turbulence to be considered a good approximation to 
spatially homogeneous decaying turbulence two criteria must be met, (i) the eddy turn- 
over time L u /u' must be small compared to the time- scale associated with the velocity 
fluctuation decay rate (du'/dx)' 1 (see also Sec. 3.3 of Townsendl 19561 ) and (ii) the rate 



of change of the turb ulent length-sc ales must be small compared to the length-scales 
themselves. Following I Corrsinl ( 1963 ) we measure 



L u du 2 L u d\ 8L U 
u 2 dx A dx : dx 

and confirm that these quantities are small for the entire decay region assessed here, i.e 
x/x* > 0.6 (figure^) and comparable with those obta ined for a regul ar grid (figure [2b) • 
Note that the 'rule-of-thumb' x/M > 40 suggested by I Corrsinl ( 1963 ) was based on the 



streamwise location where his regular grid data yielded these dimensionlcss quantities to 
be below 4%, so for our regular grid data this 'rule-of-thumb' translates to x/M > 25 and 
for our fractal square grid data to x/x* > 0.7. 

A thorough assessment of the inhomogeneity of the flow can be made by using the 
statistical equations and measuring the terms that should be zero in a statistically homo- 
geneous flow field. Starting with single-point statistics, e.g. the turbulent kinetic energy 
equation (here U\ = U, U2 = V & U3 = W denote mean flow speeds, ui = u, u 2 = v, u 3 = w 
& p are zero mean fluctuating velocities and pressure, and x\ = x, X2 = y & £3 = z are the 
components of a coordinate system aligned with the respective velocity components), 



U k dq 2 dUi d I u k q 2 u k p\ v d 2 q 2 dm dui 

o = ~ u i u j T. 7; I " l I " l 7; ^ V ~^, 7^ , (3.1) 

2 dxk oxj axk \ 2 p ) 2 dx m ox m oxk oxk 

where use is made of Einstein's notation and q 2 = u 2 + v 2 + w 2 (K = \q 2 \ over-bars 
signifying averages over an infinite number of realisations (here, over time). 

The flow statistics inherit the grid symmetries, i.e. reflection symmetry around the y 
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& z axes (as well as diagonal reflection symmetry) and symmetry with respect to discrete 
90° rotations and therefore the transverse mean velocities are negligibly small, V = W rs 0, 
and the turbulent kinetic energy equation at the centreline reduces to: 



T 



Udq 2 (— 
2 dx 



f—dU „ dU\ 8 uq 2 vq 2 \ 

[u 2 + 2uv — + 2 — 

V dx dy) \dx 2 dy 2 J 

/ d up + ^ ° vp\ + v I d 2 1 2 + 2^ 2<?2 \ 
\dx p dy p I 2\ dx 2 dy 2 J 



dui dui 
dx k dx k 



(3.2) 



where V, T, II, T> v and e are the production, triple-correlation transport, pressure trans- 
port, viscous diffusion and dissipation terms respectively. 

Data from both single- and cross- wire measurements are used to estimate all the terms 
in (|3.2[) (except the pressure-velocity correlations) along the centreline in the decay re- 
gion for Uoo = 15ms -1 (see table |4j. The pressure transport is indirectly estimated from 
the balance of (|3.2|) . The last term in (|3.2p is evaluated assuming isotropy: for the single- 
wire measurements ef^ = 15^(<9it/9x); for the cross-wire measurements one can im- 
pose one less isotropy constraint and estimate e from ef™ = 3is(du/dx) + 6v(dv/dx) 



(jSchedvin. Stegen fc Gibson It should be noted that the separation between the 



cross- wires is about 1mm and is almost 10 times the Kolmogorov length-scale so caution 
should be taken interpreting the direct measurements of dissipation using the cross- wires 
as they may be underestimated. On the other hand the isotropic estimate of the dissi- 
pation using single-wire measurements is likely to be overestimated since we show that 
(dv I dx) 2 1 (dv I dx) 2 < 2. In figure [2t the mean between the single- and cross-wire dissipa- 
tion estimates is used as the normalising quantity and the error (taken as the difference 
between the two estimates) contributes to the error bar of the normalised quantities. The 
advection 1/2 Udq 2 /dx is estimated from the non-linear least-squares power law fit of q 2 
(see Sec. 13.41 for further details) and q 2 is estimated as q 2 = u 2 (l + 2v 2 /u 2 ) with u 2 from 
the single- wire data and v 2 /u 2 from the cross- wire data; for the advection as well, the 
error is taken to be the difference between the single-wire (no anisotropy correction) and 
cross- wire estimate. The ratio between advection and dissipation can be seen (figure [5J:) 
not to be unity but tending to be approximately 1.5 beyond x/x* rs 0.8; we will return 
to this issue at the end of this subsection and in Sec. 13.41 where we estimate the decay 
rate of our turbulence. 

The longitudinal production terms are calculated from the single wire data (finer 
streamwise resolution) , whereas the transverse production terms are estimated using the 
cross- wire spanwise traverse data. The latter contribution uvdU/dy is approximately 
zero at the centreline (due to the reflexion symmetry), so it is preferred to estimate it 
just off the centreline around y ± 10mm rj L u /5 R* L v /2, to infer on its contribution in 
this region of the flow. The total contribution from the production terms around the 
centreline can b e seen (figure Eb) t o be le ss than 10% of the estimated dissipation (in 
agreement with Seoud fc Vassilicosl ( 2007 )) and beyond x > x* they become negligible 



(there is a residual production of 2-4% of the dissipation due to non-vanishing stream- 
wise mean velocity gradients). The viscous diffusion, as expected, is always negligibly 
small (table SJ. The longitudinal triple-correlation transport (table |4j shows a trend not 
dissimilar to that of the production terms, it is less than 10% closer to the kinetic energy 
peak (x/x* < 0.8) and becomes vanishingly small beyond x > x*. 

The transverse triple-correlation transport was assessed by measuring the triple cor- 
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Table 4. Turbulence statistics for five streamwise locations recorded at Uoo = 15ms~ 
Position(mm) 

x/x t 
Rex 

\fv? [ms^ 1 ] 
L u [mm] 
A [mm] 
rj [mm] 

( i+2 i) 



2 dx 

4Z = v&{£) 



~^7T 

OX 

oy 

a 

dx 2 
n d vq 2 

A a y 2 



,9V 



,a 2 ^ 

~~dy r 



(xlO 5 ) 



(xlO 3 ) 



mV 3 ] 

mV 3 ] 
mV 3 ] 



mV 3 ] 
mV 3 ] 
mV 3 ] 
\J u 2 /v 2 

(dv/dx) 2 
(du/ dx) 2 

L u f L v 



1850 


2450 


3050 


3650 


4250 


0.63 


0.83 


1.04 


1.24 


1.44 


352 


292 


253 


226 


210 


1.28 


0.99 


0.79 


0.65 


0.56 


45.7 


47.6 


50.0 


50.7 


53.6 


4.1 


4.4 


4.8 


5.2 


5.0 


0.11 


0.13 


0.15 


0.18 


0.2 


28.6 


13.9 


7.7 


4.7 


3.0 


21.6 


11.1 


6.0 


3.5 


2.2 


16.6 


8.7 


4.7 


3.1 


1.7 


0.62 


0.27 


0.11 


0.07 


0.05 


1.32 


0.12 


0.04 


0.004 


0.0007 


0.69 


0.39 


0.06 


-0.01 


-0.004 


8.22 


5.53 


3.25 


1.83 


1.08 


4.23 


1.72 


0.81 


0.42 


0.024 


3.24 


1.55 


1.29 


1.08 


0.45 


1.15 


1.13 


1.11 


1.13 


1.10 


1.39 


1.40 


1.42 


1.44 


1.46 


3.7 


3.2 


3.0 


2.7 


3.0 



relation vq 2 /2 (figure [3^) along the vertical symmetry plane of the grid (z = 0) for 
the five streamwise downstream locations specified in table [4] The transverse measure- 
ments ranged from the lower to the upper largest bars of the fractal grid (-120mm < 
y < 120 mm) and were recorded with a spacing of 20mm. The total transverse triple- 
correlation transport dvq 2 /dy (i.e. twice the transport at each transverse direction y and 
z) decreases together with the dissipation and not faster as the other measured inho- 
mogeneity terms (figure ■ It typically amounts to 40-60% of the dissipation (at the 
centreline) and perhaps surprisingly, it stays nearly the same fraction for all the assessed 
decay region. This seems to be the case not only along the centreline but for all the 
transverse measurement locations as well (figure |3Jd), although the ratio between the 
transport and dissipation arc different for different y locations and can if fact be zero 
and negative (at y/L ~ 0.35 and beyond that respectively) . 

In Sec. 13.41 wc argue that this persistent spanwise energy transport has no significant 
effect on the power law exponent of the turbulence energy decay because the dissipation 
and the lateral transport remain roughly proportional throughout the part of the decay 
region explored here. 

3.2.2. Isotropy 

The simplest assessment of large-scale anisotropy is achieved by comparing the ra- 
tio of streamwise and transverse r.m.s. velocity components, sometimes referred to as 
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-100 > > > ' > > > > 1 

0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

x/x* 

Figure 2. Homogeneity assessment for the fractal square grid in the decay region around the 
centerline for [/«> = 15ms _1 . Top figures: dimensionless streamwise inhomogeneity measures (E) 
=f ( " ) "¥§1' ( K ) llf for the ( a ) fracta l square grid (SFG), (b) regular grid (RG). Bottom 
figure: (c) T.K.E budget ()3.2|l normalised by the dissipation for the SFG at the centreline, (E) 

£ls - advection, (O) 7? . production, ( < 3) T - triple-correlation transport, (•) II - pressure 
transport, (*) £>„ - viscous diffusion. 



isotropy factor. The results of such mea surements at the c e ntrelin e are presented in 
table S] and show a fair agreement with iHurst k, Vassilicos ( 2007 ) for the same set- 
up, confirming that the flow is reasonably isotropic for all the assessed decay region, 
u'/v f « 1.1 - 1.25. T he range of isotropy factors e ncountered in our flow are comparable 
to those obtained bv lMvdlarski fc Warhaftl ( 19961) for their active grids, although further 
research shows it is possible to tun e the active grid to decrease the anisotropy of the flow 
(|Kang. Chester fc Meneveaull2003l) . Similarly it should be possible to further optimise the 
design of the fractal grids to increa se isotropy, e.g. by increas i ng the thickness rat i o as is 
suggested by the data presented bv lHurst fc Vassilicos! ( 2007 ). Hurst fc Vassilicosl ( 2007 ) 
also reported the ratio between the longitudinal and transversal integral length-scales 
{L u and L v ) for the same low-blockage space-filling fractal square grid to be L u /L v « 2, 
but this is not confirmed by the present data where the integral scales ratio is larger than 
2 as shown in table 01 even though this ratio decreases further downstream. This discrep- 
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FIGURE 3. Transverse profiles of (a) (b) at five streamwise downstream positions 

and Uoo = 15ms" 1 : (•) i/i 4 = 0.63, (□) x/x* = 0.83, (*) x/x, = 1.04, (>) = 1.24, (O) 
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Figure 4. Longitudinal and transversal one-dimensional (a) correlation function (b) energy 
spectra for x/x* = 0.63 (and for x/x* = 1.4 in the insert); (c) coherence spectra in the coordinate 
system rotated 45° with respect to the flow direction. Uoo = 15m s~ 1 . 



ancy is likely due to the calculation method of the transversal integral scales; integrating 
the transverse correlation function to the first zero crossing as lHurst fc Vassilicosl ( 20071 ) 
(incorrectly) did we recover an integral scale ratio closer to 2. 
A complementary assessment of isotropy is obtained by computing the longitudi- 
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nal and transversal correlation functions, f(r,x) = u(x)u(x + r)/u(x) 2 and g{r,x) = 
v(x)v(x + r)/v(x) 2 , and comparing g{r,x) with g tso (r) = ^d[r 2 /(r)]/dr which is the 
relation between the two correlation functions in the presence of isotropy. The compari- 
son is shown in figure 2^, for two downstream locations and it can be seen that there is a 
modest agreement between the measured and computed transverse correlation functions, 
although the agreement improves downstream. A similar comparison in spectral space is 
shown in figure |3Jd, where the isotropic relation between the longitudinal and transver- 
sal one-dimensional spectra is: F^" = [-Fii + kidFu/dki] /2. There is a fair agreement 
between the measured and computed transverse one-dimensional spectra in the 'incrtial 
region', but not at the low wave- numbers (which is consistent with L u /L v > 2) nor at high 
wave-numbers (reflec ting that (dv/dx) 2 / ( du/dx ) 2 < 2). This lack of small-scal e isotr opy 
was not reported by ISeoud fe Vassilicosl ( 20071) nor by iMvdlarski fc Warhaftl 1 19961) in 
their active-grid experiments because they filtered out the highest frequencies where their 
cross wire measurements could not be trusted. Note that in agreement with the latter 
experiments the coherence spectra (figure Hk) show that the anisotropy (inferred by the 
cross-correlation of the velocity components in a coordinate system rotated by 45°) is 
mostly contained in the large scales. The cause for this, perhaps apparent, small-scale 
anisotropy in figurc[4]and in our values of (dv/dx) 2 /(du/dx) 2 in table|4]is most probably 
the separation between the cross-wires (« Imm) being up to ten times the Kolmogorov 
length-scale. It sho uld be noted that, precis ely because of this problem, the velocity 
derivative ratios in ISeoud fe Vassilicosl ( 2007 ) were obtained for a low-pass filtered ve- 
locity signal at ki] « 0.1. This way, these authors obtained (dv/dx) 2 /(du/dx) 2 ^ 2 even 
though strictly speaking (du-i/dx) 2 = / °° k\Fu{kx) dki, where contributions coming from 
kr) > 10 _1 cannot necessarily be written off as negligible. 



3.2.3. Wind-tunnel confinement 

A qualitative assessment of the effect of flow confinement in wind-tunnel experiments 
can be made by comparing the tunnel's height/width with the flow's integral scale and 
comparing the ratio with similar experiments and with DNS. For simplicity we take the 
longitudinal integral scal^3 at the centreline to be the representative scale for each trans- 
verse section and it is typically 8.5 to 10 times smaller than the wind-tunnel width. This is 
j ust about in-line with what is typically used in DNS of deca ying homogeneous turbulence 
( Ishida. Davidson fe Kanedal Eooi IWang fe George! 120021 ). considering the boundary- 



layers on the wind-tunnels walls which reduce the effective transverse size of the tunnel 
down to 8 times the integral scale (based on the displaceme nt thickness of the boundary- 
layers) very far downstream. The active-grid experiments by Mvdlarski fc Warhaft ( 19961) 



were performed at equivalent Re\ in a similar sized wind-tunnel and produced larger 
integral-scale^ but were in line with typical decay properties and did not observe any 
of the outstanding features of our flow reported in the Subsections 13. 3[ 13.41 and 13.51 be- 
low. It is therefore unlikely that our results, namely the abnormally high decay exponent 
and the proportionality between the integral and the Taylor micro-scale, may be due 
to confinement. However it is conceivable that the effective choking of the tunnel by 
the growing boundary layers very far downstream does have some effect on the larger 
turbulence scales at these very far distances (see figure [5]). 

f For an isotropic flow the longitudinal integral-scale and the one obtained using the 3D 
energy spe ctrum ( L = n/u 2 E ( k,t)/k dk) coincide 

X Note [M vdlarski fc Warhaftl (|l996i ) used a different definition of integral-scale, but 
iGamard fc George! (|200ofi " used the same data to extract the integral-scale as defined here. 
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Figure 5. Sketch of the wind-tunnel where the decay of turbulence generated by regular and 
fractal square grids was measured. This wind-tunnel is a modified version, with an ex t ended 
test section, of the wind-t u nnel used by iHurst fc Vassilicosl (|2007l ). fSeoud fc^V assilicos (2007) 
and iMazellier k, Vassilicosl (|2010l ) in their experimental investigations of fractal generated tur- 
bulence. The boundary layers developing at the wall were estimated to have a displacement 
thickness of Si » 4mm at x = 2m (x/x* = 0.7), Si ~ 8mm at x = 3.5m (x/x* = 1.2) and 
8i » 10mm at x = 4.5m (x/x* = 1.5). 



3.3. Normalised energy dissipation rate 

It follows from this paper's introduction that for fully-developed turbulence generated 
by at least some space-filling low blockage fractal square grids, the normalised energy 
dissipation rate C £ depends both on an initial conditions/global Reynolds number Reo 
(e.g. Reo = UooX*/v) and on a local Reynolds number (Re\(x)). This distinction between 
two different Reynolds number dependencies follows from equations (1.5) and (1.7) and 
does not need to be made in the context of the Richardson-Kolmogorov phenomenology 
where the functions A and B are identical and the exponents a and (3 are both equal to 
1/2. 

The present measurements of the normalised energy dissipation rate C e for different 
Reo (by varying LT^) at two fixed streamwise downstream positions from the fractal grid 
(figure [S]) suggest that C e (Re ) is roughly constant beyond Re\(Re ) ~ 200 (figure [5]) . 
From (1.7), this observation implies that, at high enough values of Re , a = (3 and 

Ce = 15A(^-^)/B(^-^) (3.3) 

irrespetive of Reo- The facts that A is a slow- varying whereas B is fast varying function of 
^p 2 - is reflected in the steep increase of C e with x (see figure 7a). This is fundamentally 
different from the cornerstone assumption that C e is constant, an assumption which is 
approximately verified by the turbulence generated by our regular grid provided Reo is 
large enough (see figure 7a). 

The high Reo behaviour of C £ (Reo) is very comparable to that found with regu- 
lar and active- grids at similar Reynolds numbers (figure [6]) and more generally with 



other boundary-free turbulent flows such as various wakes (see e.g. lBurattini et aZ.ll2005 



Pearson et all 120021) and DNS of forced stationary hom ogeneous turbulence (see data 



compilations bv lSreenivasanlll998t iBurattini et aZ.ll2005l) . However, the fundamental dif- 



ference with the present fractal square grid-generated turbulence is that the C £ asymp- 
tote for high Reo is different for different streamwise downstream locations. This is high 
Reynolds number non- Richardson-Kolmogorov behaviour 

The key departure behind the present fractal square grid-generated turbulence be- 
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Figure 6. Normalised energy dissipation rate C £ versus Rex with Rex changing as a function 
of the inlet Reynolds number Reo for a fixed streamwise downstream location for fractal square 
grid-, active grid- and regular grid-generated turbulence. For the fractal square grid data the inlet 
Reynolds number is changed by varying the free-stream speed between 5ms" 1 < Uoo < 17.5ms" 1 
and is measured with a l w = 0.45mm sensing length single-wire at two streamwise downstream 
p ositions: (®) x/x* = 0.63 and (^) x/x* = 1.04. (♦) Active grid data is taken from table 1 
oflGamard fc George (2000) which is based on the experimental data by Mydlarski & Warhaft 
(| 19961 ) <|Gamard fc Georgg l |2000l ) computed the longitudinal and the transverse integral scales 
from the spectra, but their latter estimate yielded less scatter, hence we assume isotropy and 
use twice the tran sverse integral scale). ( <>) Regular grid data f r om th e data compilation by 
Srcenivasan (1984), figure 1 (only data bv lKistler fc Vrebalovicbl (1 19661 ) is used since no other 
experiment with more than one data point had Rex > 100). 




x(m) x(m) 

Figure 7. Normalised energy dissipation rate C £ and Reynolds number Re\ versus streamwise 
downstream location x for both the square fractal grid (SFG) and regular grid (RG) data 
recorded at different inlet velocities: (□) SFG at Uoo = 10ms" 1 , (*) SFG at Uoo = 15ms" 1 , (<1) 
RG at Uoo = 10ms" 1 , (A) RG at Uoo = 15ms" 1 , (£>) RG at Uoo = 207ns" 1 . Since the square fractal 
grid data is acquired with three different wire resolutions (see Sec. 12. 1[) we plot the arithmetic 
mean plus error bars. 
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haviour lies in the difference between the streamwise dependencies of L u /X and Rex 
(A(x) + B(x), see equations 1.4 and 1.5). For steady initial conditions (fixed Reo) 
there is a significant Re\ decrease during decay (figure EJd), whereas L u /X stays ap- 
proximately constant (figure [5^,) , leading to a steep monotonic downstream increase in 
the normalised dissipation rate C e (figure [7^,) which follows approximately the form 
C e oc (L u /X)/Rex ~ -Re^ 1 (figure [5b). Note, in particular, how the C e versus Re\ curve 
shifts to the right as Reg increases, which is clear evidence of the two independent de- 
pendencies that C e has on Re\ and Reo in this fractal-generated turbulence. 

Data were taken with probes of different spatial resolutions to confirm that these 
results are not meaningfully biased by the resolution of the measurements yielding figure 
[5J (see Sec. 12.11 for details). Nonetheless it can be seen that the lesser resolution probe 
(d w « 5/im, l w « 1mm) has a slightly lower L u /X ratio due to the underestimation of 
(du/dx) 2 , but it does not change the main observation that L u /X is effectively roughly 
constant, at least compared to the wide variation of Re\, during decay. 

We now contrast the behaviour of our fractal square grid-generated turbulence be- 
haviour with that of turbulence generated by regular-grid. Such turbulence is thought to 
follow Richardson-Kolmogorov phenomenology, although it's usually difficult to exceed 
Reynolds numbers beyond Re\ « 150 in ty pically sized la boratory wind-tunnels (at least 
if Corrsin's restriction x/M > 30 is applied, Corrsin ( 19631) ) and therefore the regular grid 
experiments are commonly at the lower end of the range of validity of the Richardson- 
Kolmogorov phenomenology. Nevertheless, our regular grid data for Uoo = 20ms~ 1 appear 
to have sufficiently high Reynolds numbers to support C s « const (figure Eb) and related 
L u /X oc Re\ (figure [8ja,) . Furthermore it can be seen that the Rex dependence of C s falls 
on the same curve regardless of how Rex is varied, whether by varying Reo or by varying 
the streamwise position of the measurement. The same observation can be made for the 
curve L u /X versus Rex- This is well-defined Richardson-Kolmogorov behaviour where 
A(x) = B(x), a = (3 = 1/2 and consequently no distinction between local and global 
Reynolds number exists. Below Rex ~ 120 direct dissipation becomes noticeable and 
causes a departure from C £ « co nst, presumably due to an insufficiently large separation 
between outer and inner scales (|Dimotakisll2000h . 

Summarising, the present fractal square grid-generated decaying turbulence is fun- 
damentally and qualitatively different from regular grid-generated decaying turbulence. 
The C e ~ const behaviour is not observed in figure |8j for the fractal square grid despite 
the moderately large turbulent Reynolds numbers Rex (around three times the Rex nec- 
essary for the regular grid to exhibit C £ = const on this plot) and the evidence that the 
global/inlet Reynolds number Reo is sufficiently large for C e to be independent of Reo 
(figure [6]). In fact the normalised dissipation rate is closer to C E ~ Rel 1 and L u /X r* const, 
which is in line with the previous experiments by lMazellier fc Vassilicosl ( 2010h . although 
the larger length of the present wind-tunnel brings to evidence that L u /X and C e Rex are 
not exactly constant in this tunnel, but are only roughly so for all the assessed decay 
region. This might be an effect brought about, perhaps paradoxically, by the eventual 
low (though not too low) values of Rex far downstream. Or it might be due to a decrease 
in the growth of L u because of the boundary layers at the tunnel walls which begin to 
have a significant thickness very far downstream in this longer wind-tunnel. As this wall 
effect might not affect the growth of A, L u /X would monotonically decrease downstream. 
Nevertheless, as we show further down in this paper, the downstream evolutions of L u /X 
and C E are consistent with a self-preserving evolution of energy spectra which can be 
made to collapse with a single-scale reasonably well, as opposed to the two different 
inner and outer scales required by Richardson-Kolmogorov phenomenology 
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Figure 8. Local Reynolds number dependence Rex(x) of (a) Integral length scale to Taylor 
micro-scale ratio L u /X and (b) normalised energy dissipation rate C e , for both the square fractal 
grid (SFG) and regular grid (RG) data recorded at different inlet velocities with single-wire 
(SW): (■) SFG recorded at [/«, = fOms -1 with a l w = 1mm sensing length SW, (□) SFG at 
Uoo = 10ms" 1 , l w = 0.45mm, (•) SFG at Uoo = 15ms" 1 , l w = 1mm, (O) SFG at t/«, = 15ms" 1 , 
l w = 0.45mm, (*) SFG at Uoo = 15ms" 1 , l w = 0.2mm, (<1) RG at Uoo = 10ms" 1 , l w = 0.45mm, 
( A ) RG at Uoo = 15ms" 1 , l w = 0.45mm, (>) RG at Uoo = 20ms" 1 , l w = 0.45mm. The dashed-dot 
lines follow the form <x Re^ 1 - The insert of the second figure is a zoomed plot of the RG data. 
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Figure 9. Uncertainty & bias due to large and small scale (an)isotropy in the observed nor- 
malised energy dissipation rate C £ behaviour. In this specific plot we re-define u and e to be 
it' 2 = u 2 (l + 2v 2 /u 2 ), e = v (du/dx) 2 [3 + 6(dv/dx) 2 /(du/dx) 2 ] and take(^) the isotropic esti- 
mates v 2 /u 2 = 1 and (dv/dx) 2 /(du/dx) 2 = 2 and (O) the anisotropy estimates of v 2 /u 2 and 
(dv/dx) 2 /(du/dx) 2 from table [4] The dashed-dot line follows oc Re^ 1 



Finally note that the large and small scale anisotropy 
u'/v' 



characterised by the ratios 
1.5, see Sec. 13.2.21 and table S} change the exact 



1.1 and (dv/dx) 2 /(du/dx) 2 
numerical values of C e and Re\ for each measurement location (see figure ^ and can be 
considered a source of uncertainty. Nevertheless, the main difference is an offset of the 
C £ versus Re\ curve and there is no meaningful change of its functional form. 



3.4. Energy decay 

The functional form of the turbulent kinetic energy decay is usually assumed to follow 
a power-law, which is mostly in agreement with the large database of laboratory and 
numerical experiments for both grid-generated turbulence and boundary-free turbulent 
flows 



(x-x y 



(3.4) 



w here v? = u' 2 . 

Mazellier &: Vassilicosl (|201dt ) proposed a convenient alternative functional form for 



the kinetic energy decay (and for the evolution of A when U^-^u 2 oc uu 2 /\ 2 is a good 
approximation) that is both consistent with the power-law decay and the exponential 
decay law proposed bv lGeorge fc Wand (l2009t l: 



2u' 2 



1 + 



1 + 



4i/a|c| 
l 2 (x )U a 

Ai/a\c\ 
P(x )U o 



-(x-x' ) 



■(x-x' ) 



n (l-c)/2c 



(3.5) 



where c < 0. In the limit of c -* it asymptotes to an exponential decay with constant 
length-scales throughout the decay, but otherwise it is a power-law decay where xq is not 
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the conventional virtual origin where the kinetic energy is singular. The two equations 
(|3.4|) & p.5j) are equivalent with n = (c - l)/2c and Xo = x' - I 2 J7oo/(4i/ac). 

Determining the decay exponent directly from (|3.4p is difficult, although feasible, since 
a non-linear fit is generally needed to determine n and x' simultaneously. For homo- 
geneous (isotropic) turbulent decaying flow where advection balances dissipation it is 
possible to obtain a linear equation for the Taylor micro-scale that can be used to de- 
termine the virtual origin, thus simplifying the task of determining the decay exponent. 
Using A 2 = I5isu 2 /e in conjunction with the advection dissipation balance characteristic 
of homogeneous isotropic turbulence (3/2 Udu 2 /dx = -e) and assuming power-law energy 
decay (|3.4[) we get 



a 10 f. 



X=—(x-x ). (3.6) 
n U 

Note that for A 2 to be linear the mean velocity has to be constant otherwise the linear 
relation holds for U A 2 . Even though advection does not balance dissipation in our fractal 
grid-generated decaying turbulence because of the significant presence of transverse en- 
ergy transport as shown in Sec. l3.2.Tl transverse energy transport and dissipation remain 
approximately proportional to each other throughout the assessed decay region and for 
the range of values of Uoo tried here. This suggests that 

f/_1^2 x (3.7) 

dx 

might be a good approximation for the decay region of our fractal-generated turbulence 
as is indeed supported by our data which show that UA 2 grows linearly with downstream 
location and even that UX 2 versus x collapses the data well for different inlet velocities 
Uoo (see figure HUh ) . 

The decay exponents of (|3.4[) and (|3.5[) are estimated using four alternative methods: 
• Method I: linear fit to UX 2 ()3.6|) to determine the virtual origin followed by a linear 

fit to t he logarithm of (13.41) to determine the exponent n, as done bv lHurst k, Vassilicos! 

(l2007h . lAntonia et al\ (|2003l ) determined the virtual origin in a similar fashion by plotting 



X 2 /(x-xq) for different xq and choosing the virtual origin yielding the broadest plateau 
(which for their regular grid experiment was Xo w 0). 

• Method II: the linearised logarithm method proposed in lMazellier fc Vassilico s (2010) 
to determine the unknowns in (13.51). 



• Method III: direct application of a non-linear least-squares regression algorithm 
('NLINFIT' routine in MATLAB™) to determine th e decay exponent and virtual origi n 
simultaneously. This is related to the method used bv lLavoie. Dienidi fc Antonia ( 2007 ). 



but allowing the virtual origin to be determined by the algorithm as well. This method 
can be applied to ()3.4[) as well as to ()3.5[) . Note that if applied to Q3.4j) as we do here, 
this fitting method does not necessarily yield a virtual origin compatible with Q3.6p . 

• Method IV: assume the virtual origin coincides with the grid location and linearly fit 
the logarithm of (|3.4[) . This crude method typically yields biased estimates of the decay 
exponent, since there is no a priori reason for the virtual origin to be zero. Nevertheless 
this is a robust method typically used to get first o rder estimates of power law decay 
exponents in many flows (e.g. the active-grid data bv lMvdlarski fc Warhaft 19961 ). 



A main difference between these methods is the way of determining the virtual origin, 
which has an important influence on the decay exponent extracted. This inherent diffi- 
culty in accurately determining the decay exponent is widely recognised in the literature 
(see e.g. iMohamed fc LaRudll990h . 
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Figure 10. Decay of turbulence generated by the regular (RG) and the fractal square (SFG) 
grid: (a) linear growth of UX 2 (b) power-law fit using method I, (c) power-law fit using method 
III (d) power-law fit using method IV. (■) SFG at U<* = 10ms" 1 , (O) SFG at [/«, = 15ms" 1 , 
( <1 ) RG at Uoo = 10ms" 1 , ( A ) RG at = 15ms" 1 (>) i RG at U«, = 20ms" 1 , (*) data from 
the Active-grid experiment bv lMvdlarski &, War haft (1993). 
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1.36 


RG 
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1.34 


0.08 
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1.25 


0.52 


1.36 


RG 


20 


1.32 


0.06 


5.47 


1.21 


0.63 


1.33 


SFG 


10 


2.57 


-0.31 


7.10 


2.51 


-0.28 


1.93 


SFG 


15 


2.53 


-0.28 


8.01 


2.41 


-0.22 


1.95 



Table 5. Decay exponents and virtual origin estimation using different methods 



The decay data for the regular grid- and fractal square grid-generated turbulence are 
well approximated by the curve fits obtained from methods I & III (sec figures [TUb & 
flOh ) and the numerical values of the exponents change only marginally (see tabled]). On 
the other hand method IV also seems to fit the data reasonably well (see figures ITUa) 
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Figure 11. Turbulent kinetic energy decay of turbulence generated by the fractal square grid 
fitted to (|3.5[1 using method II (dashed-dot line) and method III (solid line). The data range 
used in method II is 0.6 < x/x* < 1.1, which corresponds to the streamwise region assessed 

by Mazcllicr & Vassilicos (2010]). Streamwise data was taken at two fixed inlet velocities: (■) 
[Too = 10ms -1 , (O) U<x> = 15ms -1 . Notice that for 0.6 < x/x* < 1.1 the two methods appear to fit 
the data reasonably well, but further downstream the differences become evident. 



but the exponents retrieved for the fractal grid data are n » 2, slightly lower than the 
exponents predicted by the other methods n « 2.5. The virtual origin which is forced 
to xo = in method IV leads to a slight curvature in the log(u 2 ) versus log(x) data 
(almost imperceptible to the eye, compare the fractal grid data in figures HUb & HUH) 
and a non-negligible bias in the estimated exponents. Nevertheless the difference in the 
power laws describing the measured regular grid- and square fractal grid-generated tur- 
bu lence is quite clear. Fo r comp leteness, the results from the experimental investigation 
by iMvdlarski k, Warhafti ( 19961 ) on decaying active grid-generated turbulence are added 



in figure HUH . They applied a fitting method eq uivalent to method IV and reported a 
power-law fit yielding a decay exponent n = 1.21. Kang et al. ( 20031 ) employed the same 
method to their active grid-generated turbulence data and retrieved a similar result, 
n= 1.25. 

Note that there are residual longitudinal mean velocity gradients (which cause a resid- 
ual turbulence production of about 3% of the dissipation, see Sec I3.2.ip and therefore it 
is preferred to fit u 2 data rather than u 2 /U 2 data. Nevertheless we checked that fitting 
u 2 /U 2 data does not meaningfully change the results nor the conclusions. 

Concerning method II it can be seen (table [5]) to be the most discrepant of the 
four methods yielding a mu ch larger decay exponent. This method was proposed by 
Mazellier fc Vassilicos! ( 201 d( ) to fit the general decay law Q3.5P and is based on the lin- 
earisation of the logarithm appearing in the logarithmic form of (|3.5p . i.e. 



log( 



U ) 



log 



l + c 
~2c 



log 1 + 



Avac 



<*))■ 



(3.8) 



Linearisation of the second logarithm on the right hand side of 



assumes 



4mc / 



x'q) « 1. This quantity, as we have confirmed in our data, is indeed smaller than unity 
and for the farthest position 4vac/(l 2 (xo)Uoo)(x - x ) w 0.3, but the fact that this lin- 
earised method does not yield results comparable to methods I and III suggests that 
the linearisation of the logarithm may be an oversimplification. In figure [TT] the kinetic 
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energy decay data of turbulence generated by the fractal square grid is shown along with 
the fitted curves obtained from methods II and III in a plot with a logarithmic ordinate 
and a linear abscissa. In figure [TTa the data taken at pos itions beyond x/x* ~ 1 . 05 are 
excluded in order to compare with the results presented in iMazellier fc Vassilicosl (j2010h 



where the data range was limited to 0.5 < x/x* < 1.05. Visually, in figure [TTk. the two 
different fitting methods appear to fit the data reasonably well and thus the linearisation 
of the logarithm in (|3.8p is justified in this limited range. Note, however, that the two 
fitting methods yield very different decay exponents because they also effectively yield 
different virtual origins: for example at Uoo = 15ms -1 method III yields (l + c)/(2c) « -2.4 
whereas method II yields (1 + c)/(2c) « -8.0. In figure [TTb. where no data is excluded, 
it can clearly be seen that the two methods produce very different curves and very dif- 
ferent decay exponents (note however that the use of a longer test section, which allows 
the assessment of the decay behaviour further downstream, comes at the cost of having 
thicker boundary layers developing at the walls which can have an increasing influence 
on the largest turbulent eddies, as discussed in Sec. I3.2.3|) . 

3.4.1. Influence of transverse transport on power-law decay exponent 

It is shown in Sec. 13.2.11 that dissipation does not balance the advection but that 
the two are roughly proportional throughout the measured decay region of the fractal 
square grid-generated turbulence. It is also shown in that section that this imbalance is 
mostly due to transverse triple-correlation transport which remains roughly 50 - 60% of 
the dissipation throughout the measured region (with no clear increasing or decreasing 
trend) , whereas turbulence production and longitudinal triple-correlation transport terms 
become negligible well before x < x*. Pressure transport, calculated from the kinetic 
energy balance, may also play a noticeable role of countering a fraction (typically between 
1/4 and 1/3) of the triple-correlation transport. Based on these results, equation (|3.2[) 
which holds at the centreline reduces to 



Udq 2 

= -e + 

2 dx 



dy 2 



(3.9) 



The decay rate of the kinetic energy as the turbulence is advected downstream (effectively 
the advection term) is now determined both by viscous dissipation and by a net effect of 
removing energy from the centreline and transporting it to the sides. As in the portion 
of the decay region of the fractal-generated turbulence where we take measurements this 
loss rate to the sides remains approximately proportional to the dissipation rate, i.e. 

Udq* 

2 dx ~ X£ 

where x = 1 + [ dvq 2 dy -IL]/e * 1.5 (figure [2b), we can expect the decay exponent n to be 
set by the dissipation rate e (irrespective of what sets the dissipation rate). Indeed, the 
higher power law decay exponents exhibited by the fractal-generated turbulence can be 
accounted for by the fact that C e ~ Re^ 1 (see Sec. l3.3p and consequently the steep increase 
of C e = eL u /u' 3 with strcamwise location. In other words, an increasing proportion of 
u' 3 /L u is being dissipated at increasing streamwise locations which leads to an increase 
in the power law decay exponent relative to the C E = const case. 

In figure [T2] we plot in logarithmic axes the streamwise decay of the advection, the 
dissipation and the transverse triple-correlation transport (which are all measured inde- 
pendently) and they indeed seem to follow straight lines (i.e. power laws) with the same 
slope (i.e. power law exponent), thus supporting our argument. 

To further substantiate our argumentation one more set of experiments were con- 
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Figure 12. Downstream decay of: (•) kinetic energy ¥r^§—, ('-') dissipation e, (*) transverse 



triple-correlation transport 2 



d_ I vq 2 \ 
dy \ 2 )> 



for U = 15ms 



ducted. Anemometry measurements at an inlet velocity of Uoo = 15ms -1 using a l w = 
0.5mm sensing length single-wire were recorded between 0.63 < xjx* < 1.44 along four 
parallel lines aligned with the mean flow and crossing the grid at z = 0, 0, 20, 20 mm 
and y = 0, 40, 80, 120mm (z = is the vertical plane of symmetry of the grid). From the 
transverse triple-correlation transport measurements for 2 = (figure [3Jd) we expect the 
contribution from this term to be very different at the centreline (where it is maximal) 
and off the centreline where is can be roughly zero (y « 80) or negative (y > 100). How- 
ever, if a value of x can be defined that is constant throughout the streamwise decay 
range assessed here for each transverse (y, z) position, then the argument outlined in the 
previous two paragraphs will hold even if \ varies with transverse positions, as indeed it 
does. The consequence is that, in the decay region assessed, the decay exponent n should 
remain about the same at all these transverse positions and also remain unusually large 
due to the C e behaviour. The data for the different transverse locations are fitted using 
method III and the results (see table [5]) are encouraging. In spite of some variation in 
the best fit power-law decay exponents, the numerical values of these exponents are all 
relatively close to each other ranging between 2.3 and 2.6. We note that these exponents 
are larger than all boundary- free turbulent flows listed in table [TJ 

Finally, as some presence of turbulence production and longitudinal transport remains 
for some distance downstream of x pea k ~ 0.45a;* (though not in any significant way be- 
yond a;*) we explore how the power-law fits of the turbulence energy decay change when 
the smallest streamwise location considered in the fit is increased. We do this both for 
centreline and off-centreline data and report our results in figure 1131 On the centreline 
the decay exponent and virtual origin remain approximately the same within the scatter 
(n « 2.4, xq/x* « -0.3), but they show a respectively decreasing/increasing tendency off- 
centreline up to x/x* « 0.8. At any rate, the decay exponents n > 2.0 for all our data. 



In conclusion the decay exponents for the present fractal-generated turbulence mea- 
sured both at the centreline and off the centreline in the region 0.6 < x/x± < 1.5 are 
consistently higher than those in all boundary-free turbulent flows listed in table 1 and 
much higher (by a factor betwe en 4/3 and 2) than those of decaying turbule nce gener- 
ated by regular and active grids ( Mvdlarski fc Warhaftlll996llKang et al. 2003f l. It might 
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Table 6. Decay law estimates along four parallel streamwise oriented lines at the centreline 
and off the centreline between 0.63 < x/x+ < 1.40 obtained from method III. 



y (mm) n xo/x* 



2.42 -0.27 

40 2.61 -0.29 

80 2.27 -0.11 

120 2.63 -0.39 



3 I 1 ! ! ! ! 1 0.5 

2.8 

2.6 
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1 I i i 1 i i 1 -1 I i i 1 i i 

0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.6 0.65 0.7 0.75 0.8 0.85 0.9 

Figure 13. Decay law estimates for different data ranges x m i n < x/x t < 1.40, U = 15ms -1 : (•) 

y = 0mm, ('-') J/ = 40mm, (*) y = 80mm, y = 120mm. (a) decay exponent n and (b) virtual 
origin xo/x* obtained from method III 



be interesting to note that in many boundary-free turbulent flows a conserved quantity 
such as u' 2 L M+1 = const exists. Look at tabic Q] and note that M = 1, 3, 5, 7 for the four 
wakes, M = -1 for the mixing layer, M = 0, 1 for the jets and M ^ 2 for regular-grid 
turbulence. If the flow is also such that Udu' 2 /dx <x -e then C e = const implies 

2(M+1) 

n = 

M + 3 

and C e ~ Re^ 1 implies 

M+ 1 



(which is larger than n = 2(M + 1)/(M + 3) provided that M > 1). Considering, for 
example, the range M ^ 2, the exponent n corresponding to C e ~ -Re^ 1 is at least 5/4 
times larger than the exponent n corresponding to C e = const, and is generally much 
larger. If M = 3 or M = 4 then C e ~ Re^ 1 implies n = 2 or n = 2.5, close to what is 
observed here, whereas C £ ~ const implies n = 4/3 or n = 10/7. 

At this stage we do not have any proof that a conserved quantity such as u' 2 L AI+1 = 
const exists for our fractal-generated turbulence. The previous paragraph is therefore 
only indicative and serves to illustrate how a C e which is a decreasing function of Re\ 
can cause the decay exponent to be significantly larger than a C e which is constant dur- 
ing decay and can even return decay exponents comparable to the ones observed here. Of 
course the decaying turbulence we study in this work is not perfectly homogeneous and 
isotropic because of the presence of transverse turbulent transport of turbulent kinetic 
energy and therefore significant gradients of third-order one-point velocity correlations. 



28 



P. C. Valente and J. C. Vassilicos 



As a consequence, a conserved quantity such as u' 2 L M+1 = const, if it exists, cannot result 
from a two-poi nt equation such as the von-Karman-Howarth equation for homogeneous 



turbulence (see IVassilicosI 120111) . We leave the investigation of conserved quantities in 
third-order inhomogeneous decaying turbulence such as the present one for the future 
(we include gradients of pressure- velocity correlations in the term " third-order inhomo- 
geneous"). 

Nevertheless, it is clear that the dissipation rate of kinetic energy is increasingly larger 
than u' 3 /L as the turbulence moves further downstream in cases such as the present 
one where C e increases in approximate proportion to l/Re\ as the turbulence and Re\ 
decay. In the absence of any other type of loss or gain of kinetic energy, and assuming 
no counter-effect of C £ on the integral scale, a much steeper decay (e.g. much larger 
exponent n) will result than if C e was constant during decay. In the present case where 
loss of energy also occurs by turbulent transport, see equation ()3.9|) . this conclusion can 
remain the same in the region assessed only if, in that region, the loss of energy by 
turbulent transport remains proportional to the loss of energy by dissipation, as indeed 
observed. 

The question then naturally arises whether this balance between turbulent transport 
and dissipation persists for the entire decay range all the way to very large values of x/x*, 
much larger than those accessible here. If it does, then the implication is that perfectly 
homogeneous isotropic turbulence is impossible at any stage of the decay. If it does not 
and if turbulent transport starts to decay much faster than dissipation beyond a certain 
x/x*, then a turbulence that is third-order homogeneous and isotropic may well appear 
if it has the time to do so before the final stages of decay. If C e continues to increase 
nearly as l/Re\ in such a third-order homogeneous isotropic turbulence then the decay 
will remain exceptionally fast with values of n such as the present ones. However, it 
may be that the unusual behaviour observed here for the dissipation rate e (a two-point 
statistic) is in fact the result of gradients in particular one-point statistics such as third- 
order velocity correlations and pressure- velocity correlations, i.e. inhomogeneities. Either 
way, the consequences can be far reaching and call for much future research, in particular 
re-examinations of Reynolds number dependencies of C £ in all manner of turbulent flows, 
in particular boundary- free turbulent flows such as those listed in table [1] 

As a final remark, note that the data points for transverse turbulent transport in figure 
[12] seem to curve downwards at high x/x* . However we cannot extrapolate much from this 
observation as we do not measure pressure directly and we do not know how gradients 
of pressure- velocity correlations curve at high x/x*. 

3.5. Collapse of the energy spectra and structure functions 
As explained in Seoud fc Vassilicosl ( 2007 ) and Mazellier fc Vassilicos! ( 20 loh single- length- 



scale self-preserving energy spectra can allow for L u /X = const during decay. This can 
be assessed by plotting the normalised energy spectra for different positions along the 
mean flow direction and evaluating the collapse of the data or the lack thereof. It should 
be mentioned that the three-dimensional energy spectrum and one-dimensional energy 
spectra can be shown to be equivalent for an isotropic flow. It should also be noted that 
the flow is not exactly isotropic as discussed in Sect. l3.2.2| so we might expect some effect 
on the spectral collapse. 

3.5.1. One- dimensional energy spectra 

We begin by illustrating the qualitative difference between the collapse of the nor- 
malised energy spectra (using large scale variables: u' 2 , L u ) of turbulence generated by 
the regular grid and by the fractal square grid, see figure El The data for the regular 
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Figure 14. Compensated ID energy spectra, normalised with u 2 , L u at two streamwise locations 
for the (a) regular grid- and (b) fractal square grid-generated turbulence data. Both plots have 
roughly the same Reynolds number ratio, Re\ 1 /Rex 2 ~ 1.3 (see Appendix O. The data are 
recorded at Uoo = 20ms~ 1 and [/«, = 10ms -1 respectively. 



grid are taken in a region where L u /X oc Rex and C e » const, see figure [51 The normalised 
spectra measured in the lee of the regular grid show a good collapse at the low frequen- 
cies but not at the high frequencies, which is in-line with Kolmogorov's theory. On the 
other hand it can be seen that the normalised turbulence spectra generated by the frac- 
tal square grid appears to collapse at all frequencies, in-line with t he singlc-lcngth-scale 
assumption as previously observed bv lMazellier fc Vassilicosl (|20ld ). 

In order to complement the previous results, the assessed decay region is extended 
allowing to further test the single-length-scale assumption. The normalised spectra of 
decaying turbulence downstream of our fractal square grid are shown in figure [15] using 
both the integral-scale and the Taylor micro-scale. For the extended region it can be 
seen that the normalised spectra using the Taylor micro-scale do collapse for the entire 
frequency range, although the collapse using the integral-scale at high frequencies is 
modest for kL u > 40 where the furthermost point (x/x* = 1.41) is taken into account. 
However the discrepancy between data at Re\ = 324 and Re\ = 210 is much too small 
compared to the lack of collapse which would occur if the data obeyed Richardson- 
Kolmogorov scaling as in figure [T4k. It should be noted that in theory the collapses with 
L u or with A should be identical if L u oc A, but as was seen in figure [8J this is not verified 
exactly in our wind-tunnel's extended test section. 

Nevertheless, in Appendix lAl we propose a methodology for making a rough estimate 
of the quality of collapse of normalised spectra at high frequencies and we find that it 
depends on the logarithm of the Reynolds number ratio Re\ t l Re\ 2 at two streamwise 
distances x = £i and x = £2 with a pre-factor which depends on the behaviour of L u /X 
during decay. In the Appendix, we apply this methodology to the active grid data of 
Kang et al. ( 20031 ). for which there is evidence of a Richardson- Kolmogorov cascade, and 
show how spectral collapse with outer variables can be misleading because the Reynolds 
number ratio is small. The same methodology applied to our data shows that we are not 
fully able to conclude on the very high frequency end of fractal grid-generated energy 
spectra. 

Pre-multiplying the ID energy spectra by the square of the frequency yields the Fourier 
spectrum of (du/dx) 2 , so a second test to the single-length-scale assumption is to assess 
the collapse of this isotropic equivalent of the dissipation spectra. The data, plotted in 
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Figure 15. Compensated ID energy spectra of turbulence generated by the fractal square grid 



at three streamwise downstream locations {Re\ = 324, 273, 210) at U a 
(a) u 2 and L u (b) u 2 and A (c) v and e. 



15ms , normalised by 




Figure 16. Compensated ID energy spectra of turbulence generated by the fractal square grid 
at three streamwise downstream locations corresponding to Re\ = 324, 273, 210 at Uoo = 15ms _1 , 
normalised by (a) u 2 and L u (b) u 2 and A. 
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figure ITBI show a reasonable collapse onto a single curve using both length-scales, though 
it can be seen that the peak of the pre-multiplied spectra does not collapse perfectly. 
This may, to some extent, be an effect of the slight anisotropy of the flow, since it affects 
the large scale variables the normalisation is based on. It is shown in the following section 
how it is possible to partly account for this effect by computing the three-dimensional 
energy spectrum. 

3.5.2. Three- dimensional energy spectra 

The 3D energy spectrum is computed using the two-component velocit y signal from the 



cross- wire measurements, with a similar algorithm to the one presented in lHelland fc Van Atta 



(|1977I ). The central assumption of the algorithm is isotropy in order to relate the one- 
dimensional total energy spectrum Fu(ki) = Fn(ki) + ^22(^1) + ^33(^1) with the three- 
dimensional spectrum E(k), 

E(k) = -k^ (3.10) 
dk\ 

where the transverse one-dimensional spectra are considered to be approximately the 
same, i.e. ^22(^1) ~ -^33(^1)- The first derivative of the spectrum is computed using the 
logarithmic derivative proposed bv lUberoil (|l963h : 



am k\ 

The 3D energy spectrum was evaluated at 50 logarithmically spaced frequencies, and 
a 2 nd -order polynomial was fitted between two neighbouring frequencies using a least- 
squares-fit in order to obtain a smooth derivative of the spectrum. 

From the 3D energy spectrum the integral scale L, the turbulent kinetic energy and 
the Taylor micro-scale can be recovered. The difficulty in accurately determining the low 
frequency range of the energy spectra and consequently estimating the integral length 
scale should be noted. For this reason, the assessment of the spectrum's slope near k -> 
was not possible. 

The normalised compensated spectra are shown in figure 1171 while the normalised 
enstrophy spectra are shown in figure 1181 It is rewarding to see that the collapse of the 
3D energy spectrum presents less scatter than the ID spectrum thus offering support to 
the self-preserving single-length behaviour of turbulence generated by the fractal square 
grid. Hence, some of the deviation from single-scale self-similarity collapse of the ID 
spectra in figures [15] & [16] is due to the moderate level of anisotropy present in the 
turbulence. 

3.5.3. Second-order structure functions 

The collapse of the second order structure functions using u' 2 and L u is shown in figure 
119a . Similarly to what has already been discussed for the spectra, this structure function 
collapses well at both low and high separations in the case of our fractal-generated 
turbulence. However, this is clearly not the case for the turbulence generated by the 
regular grid (see figure [T^b). 



4. Conclusions and issues raised 



The decay of regular grid- and fractal square grid-generated turbulence have been ex- 
perimentally investigated using constant temperature hot-wire anemometry. The main 
contribution of the present work is to complement previous research on the decay of frac- 



tal grid-generated turbulence (e.g. iHurst fc Vassilicoa 120071 : ISeoud fc Vassilicosl 12007 : 
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Figure 17. 3D energy spectra of turbulence generated by the fractal square grid at three 



streamwise downstream locations corresponding to Re\ 
normalised by (a) u 2 = 2/3q 2 and L (b) u 2 and A. 
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Figure 18. Enstrophy spectra of turbulence generated by the fractal square grid at three stream- 
wise downstream locations corresponding to Re\ = 300, 238, 210 at [/«, = I5ms~ and normalised 
by (a) = 2/3^2 and L ( b ) ^2 and A 
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Figure 19. Second-order structure function at three streamwise downstream locations nor- 
malised by u 2 and L u : (a) SFG recorded at Uoo = 15ms" 1 , Re x = 323, 273, 210 for 
x/x, = 0.69, 0.90, 1.41 (b) RG recorded at Uoo = 20ms" 1 , Re x = 156, 137, 120 for x/M = 22, 38, 68. 
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Mazellier fc Vassilicod 2010h by doubling the extent of the assessed decay region with 



the aim of investigating the persistence (or lack thereof) of the reported high decay ex- 
ponents and the suppressed Richardson-Kolmogorov cascade. The present experimental 
investigation also complements the previous research by studying the effect of the hot- 
wire spatial resolution, carefully assessing the homogeneity of the flow during decay and 
taking anisotropy into account in the energy spectra. 

We find that for streamwise downstream positions beyond x/x* w 0.6 the turbulence 
is close to homogeneous except for a persistence of pressure transport and transverse 
energy transport and decays such that L u /X rj Const whilst Re\ sharply decreases, at 
least up to the furthermost downstream position investigated. However L u /X increases 
with increasing grid Reynolds number, e.g. Re$ = Uo pt^jv. This observation is i n di- 
rect conflict with the Richardson-Kolmogorov cascade ( Mazellier fe Vassfflcodl20ich . be- 
lieved to be dominant at this range of Taylor-based Reynolds numbers Re\ in various 
boun dary-free turbulent flows, including regular gr id- and active grid-generated turbu- 
lence ( Burattini et al. 2005t Sreenivasanll 1 984l 1 1 998 ) . It must be noted, however, that the 
vast majority of existing data is taken at fixed streamwise locations x and varying inlet 
Reynolds numbers Reo and as we show in section 13.31 for fractal grid-generated turbu- 
lence, the streamwise downstream Reynolds number dependence Re\(x) isn't necessarily 
the same. 

We observe that the energy spectra and the 2 nd order structure function are much 
better described in the p resent fractal square grid-generated turbulence by a single-scale 
self-sim ilar form than by iKolmogorov ( 194ll ) phenomenology. Note that by iKolmogorov 
(194l|) phenomenology we mean, not only the necessity of two dynamically relevant 
sets of variables, outer and inner, that collapse the low- and the high-frequency part of 
the spectra respectively, but also that L u /X oc Re\ and C £ = Const, which implicitly 
dictates the rate of spreading of the high-frequency part of the spectra normalised by 
outer variables and vice-versa. That tu rbulence generated by the present fractal square 
grid does not obey Kolmogorov ( 194lh phenomenology is clear, for example, from the 
comparison between figures IT4k andlMb. 

We also confirm the observations of lHurst fe Vassilicosl ( 2007 ) and Mazellier fe Vassilicos 
( 201Cf ) concerning the abnormally high power- law decay exponents, compared with most 
boundary-free turbulent flows (see table[TJ , in particular regular and active grid-generated 
turbulence {nsFG » tirg^^ag by a factor between 4/3 and 2), and we confirm their 
persistence further downstream (at least up to x/ x* 1.5). However, ou r resul ts do not 
support the view in Hurst fc Vassilicos ( 2007 ) and iMazellier fc Vassilicos ( 2010l ) that the 
turbulence decay is exponential or near-exponenti al. We infer, by comparing o ur experi- 
mental results with the active- grid experiments of iMvdlarski fc Warhaftl (jl996j), that the 
reason for the very unusual turbulence decay properties generated by the fractal square 
grids cannot be a confinement effect arising from the lateral walls. The two experimen- 
tal investigations report completely different turbulence properties during decay, even 
though both experiments were performed on a similar sized wind-tunnel and, in fact, 
the integral length-scales generated by our fractal square grid are typically less than 
half the integral length-scales generated by the active-grid. Our fractal-generated tur- 
bulence is third-order inhomogeneous in the sense discussed in subsection 13.4.11 but, to 
our knowledge, no homogeneity studies of active grid-generated turbulence exist to this 
date which are as thorough as the one presented here, and it is therefore not possible 
to fully compare homogeneity and isotropy levels of the two types of turbulence. The 
presence/absence of turbulent transport of pressure and kinetic energy have not been 
investigated in sufficient detail in either active or regular grid-generated turbulence and 
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it remains unknown to what degree and how far downstream these types of turbulence 
are third-order homogeneous and isotropic. 

Although we find a general agreement with the previous results on the decay of fractal 
grid-generated turbulence, some new issues are raised by the present experimental results 
due to the extended wind-tunnel test section. We find that L u /X is in fact not perfectly 
constant, but slowly decreases with Re\, and that the spectral collapse using large scale 
variables is not perfect at very high wavenumbers as it ought to be for exact single- 
scale self-preserving turbulence decay. Possible causes for these two observations will be 
investigated in future work and include: (i) small scale corrections to the single-scale self- 
preservation, (ii) moderately low Reynolds number limit to the validity of single-scale 
self-preservation and (iii) excessive thickness of the confining wall boundary layers far 
downstream interfering with the growth of the largest eddies of the turbulent flow due 
to insufficient ratio between the wind-tunnel width and the integral length scale. 

As a final remark we note that the study of freely decaying turbulence requires exper- 
iments where (i) a wide range of i?eo can be achieved by modifying the initial conditions 
and (ii) a wide range of Re\ values must be straddled during decay. This is emphasised 
by the analysis presented in Appendix [X] where we give quantitative criteria for truthful 
spectral collapse and where we show, in particular, that whereas active grids gener- 
ate high Reynolds numbers they also generate a narrow logarithmic range of Reynolds 
numbers during decay thus making it impossible to confirm the Richardson-Kolmogorov 
cascade via spectral collapse. 

So far, only modest ranges of Rex during decay have been achieved with regular grid- 
and active grid-generated turbulence due to the typically slow decay rates of the tur- 
bulence they generate. The fractal square grid-generated turbulence offers the unprece- 
dented possibility of generating high intensity decaying turbulence with a very wide range 
of Re\ values during decay and approximately homogeneous mean flow and turbulence 
intensity profiles. 

We are grateful to Charles Meneveau, Beat Liithi and one anonymous referee for nu- 
merous suggestions which have helped us to very significantly improve the paper. P.C.V. 
acknowledges the support of the Portuguese Foundation for Science and Technology 
(FCT) under grant number SFRH/BD/61223/2009. The support from Dantec Dynamics 
Ltd (UK) in the design of the sine-wave testing hardware is also acknowledged. 



Appendix A. A note on the energy spectra collapse in turbulence 
generated by active grids 

Active-grid experiments can generate relatively high R e\ turbulence in a typically 
sized laboratory wind-tunnel ( Mvdlarski fc Warhaft 19961 ) and are thus a potentially 
good test case to compare the fractal-generated turbulence with. The comparison of the 
energy decay rate was shown in Sec. 13.41 and here we focus on the co llapse of the energy 
spectra. The data is taken from the experimental investigation by iKang et al ( 20031) 



on the decay of active grid-generated turbulence. In their paper the longitudinal energy 
spectra for four downstream positions is provided in tabular form and we use this data 
here to create the plots in figure [20] 

The pre-multiplicd longitudinal energy spectra (figure I20[) are normalised using both 
outer (u 2 and £) and inner (e and 77) variables. Note that I is the pseudo-integral-scale 
defined as I = 0.9u' 3 /e which is proportional to the integral scale if and only if the 
dissipation coefficient C £ is constant during decay, in this case C £ = 0.9. 

At a first glance the results seem striking since both outer and inner variables seem to 
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Figure 20. Compensated ID energy spectra at four streamwise locations 
(Rex = 716, 676 650, 626) a t U_~ 11ms -1 , normalized by (a) up and £ = 0.9u' 3 /e (b) s 
and n. Data from lKang et al] (|2003f ). 
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Figure 21. Sketch of two spectra at two streamwise positions x = £i and x = £2 normalised 
with outer variables spreading at high frequency. 



be collapsing the spectra. Thus one could conjecture that turbulence generated by active 
grids is self-similar and has only one determining length-scale. This would in fact be the 
case if the Reynolds number Rex remained constant during the turbulent kinetic energy 
decay and consequently L oc A oc 77, which is the first ever self-preserving decay proposed 
( Karman fc Howartlj[l938h . ] nstead Rex decreases during decay. Hence this collapse can 
only be apparent, but not real. 

A simple method of estimating the necessary range of Reynolds numbers Rex for the 
collapse to be meaningful is now presented where it is shown that the collapse (or spread) 
of a normalised spectrum at two streamwise locations is only significant if the logarithm 
of the respective Reynolds numbers' ratio is large, typically \og(Rex 1 /Re\ 2 ) > 1/4. The 
starting point in this methodology is the assumption that a given scaling is correct (e.g. 
Kolmogorov or single-length scalings) which then allows the quantification of the spread 
for a given Rex range of any other attempted normalisation. 

We outline the method by considering the dissipation range of the longitudinal spec- 
trum and assuming the Kolmogorov scaling is correct, i.e. Fn(k,x) = e 2 ^ 3 rj 5 / 3 f(krj), but 
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this methodology is easily extendable to the energy containing range of the spectrum as 
well as to the case where the single-length-scaling is correct. 

Consider two streamwise distances x = £1 and x = £2 and write r\\ = r\2 = 

»/(&), Ai = A(£i), A 2 = A(&), £1 = £2 = £*(&), «i = «'(£), 4 = «'(&), 

£i = e(£i), £ 2 = £(£2) for the Kolmogorov scales, Taylor micro-scales, integral scales, 
r.m.s. turbulence velocities and dissipation rates at these two locations. We take £2 > £1 
so that Re[ Jl = u^Li/v > RtL 2 = u' 2 L 2 /i/. 

Choose two wavenumbers k\ and /«2 in the dissipation range such that fcir?i = k2f}2 and 
/(fci?7i) = .f(fc 2 ?7 2 ) by assumption. If one would normalise the same spectra in this range 
using u' 2 & L u , the dependence of the normalised spectra on x would explicitly resurface, 
i.e. Fn(k,x) = u' 2 L u f*(kL u ,x) (see figure[2T]). Since e = C e u l3 /L u with C E independent 
of x in the Richardson-Kolmogorov phenomenology, it follows that L u /rj = Cl^Re 3 / 4 and 
it is possible to show that 



nhL^) = r(fc 2 £ 2 ,6) = r{ - k2l2 ^ 2) (fSf) ' 



(Al) 



and 



so that f*{k\Li,^i) + /*(fc 2 L 2 , £2) and kiL\ + k 2 L 2 . 

We define the spectral spread which characterises the degree of non-collapse by the 
form Fu = u' 2 L u f*(kL u ) as 

= logCMi) - log(k 2 L 2 + 6kL 2 ), (A 3) 

where = f*{k 2 L 2 + SkL 2 ,t; 2 ), see figure [2T1 There are two contributions to 

the spectral spread, one from the rescaling of the abscissas, fciLi + k 2 L 2 , and another 
from the rescaling of the ordinates. From Rcl 1 > R&l 2 and equations (|A (|A 2|) we 
know that fci-Li > ^2-^2 ana - /*(^i£i>£i) < /*(^2-k2, £2) so that the two contributions 
to the spectral spread can actually in principle, cancel each other. However, the second 
contribution depends on the functional form of f*(kL u ) and therefore it is not possible to 
quantify its spectral spread contribution without an analytical expression for f* (kL u ,^). 
Nonetheless, as is shown below, we can estimate a bound for this contribution, so that 
in the end we can estimate a upper and lower bound for the expected spectral spread "J 
characterising the degree of non-collapse by the alternative scaling. 

The contribution to the spread \& from the abscissa's rescaling alone (which is the 
upper bound) is given by (using (|A 2ft ) 

* max = logCfclx) - log(fc 2 i 2 ) = Wfj 1 ) = 1 1o 4rJ^) ■ ( A4 ) 

(for the last equality, ([1.6)1 was used to relate the integral scale to the Taylor micro-scale 
with C e = Const from Richardson-Kolmogorov phenomenology). 

The contribution to the spread $ from the ordinate's rescaling is measured as a fraction 
of the abscissa's rescaling 

$ _ \og{k 2 L 2 + 5kL 2 ) - \og(k 2 L 2 ) 
log(fciLi) - log(fc 2 £2) 

so that $ = for 5kL-2 = (ordinate rescaling has no effect) and $ = 1 for 5kL>2 = 
k\L\ - fc2-^2 (ordinate rescaling cancels the abscissas rescaling). It is possible to show 
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using a first order Taylor expansion in logarithmic coordinates that we can re-write the 
function $ to leading order as 



aio g (/*(iog(fcL u ),o) 



d \og{kL u ) 



(A 6) 



fcL„ =k 2 L 2 



Since the spectra in the dissipation range roll-off faster than any power law we can 
always find a high enough wavenumber kt(p)L u for which the tangent of the spectrum 
(in logarithmic coordinates) is steeper than (kL)~ p given an exponent p (see figure l2Tj) . 
Consequently, for a given choice of p, we get an upper bound for $ for wavenumbers 
above k t (p)L u which is $ ma x = 5/(3p). Therefore we can estimate a lower bound for the 
spectral spread as ^ m in = ^max ~ ^max an( l thus 

3, (Re Xl \ 5 _ 3. (Re Xl \ , A 

In figure l2"0k we plot the (kL u )~ p function with p = 4 and it can be seen that for 
wavenumbers higher than ktL u » 400, the tangent of the spectra (in logarithmic coordi- 
nates) are steeper. Hence, according to (|A 7J) . for k t L u > 400 (taking into account the Re\ 
range of the data presented in figure [20k) the spectral spread will be around 5% < \P < 9% 
of a decade, which can easily be confounded with scatter. Ther efore, the apparen t col- 



lapse observed in the spectra from the active-grid experiments of lKang et al\ (|2003l) (see 
figure [20k ) may be misleading as it is simply the result of a small range of Re\ variation 
during decay (from 716 to 637). Note that the same misleading collapse occurs in the 
low wavenumber range of the spectra plotted in figure 120b where the Kolmogorov inner 
variables were used for the normalisation. 



We can repeat the exact same analysis for the case where the Richardson-Kolmogorov 
cascade is suppressed and assume the validity of Fn(k,x) = u' 2 L u f*{kL u ) and L u /X « 
Const. We can then carry out the same calculation as above to obtain the spectral 
spread = log(fcir/i) - log(fc 2 ?72 + Shift) for k\Li = k 2 L 2 when attempting to collapse the 
data with Kolmogorov variables. We would then obtain the same expression to quantify 
the spread contribution caused by the rescaling of the ordinatcs relative to the total 
spread, (|A 5[) with a suitably redefined $ where the outer scales L u have been replaced 
by inner scales r\. Note that spread of the high frequency spectra normalised by outer 
variables assuming that F n (h,x) = e 2/3 rj 5 / 3 f(kn) holds is the same as the spread of 
the high frequency spectra normalised by Kolmogorov inner variables assuming that 
Fn(k,x) = u' 2 L u f*(kL u ) holds. Hence, using L u ~ A, the spectral spread resulting from 
an attempt to collapse with Kolmogorov inner variables spectra which obey complete 
self-similarity is 



* = log(fci/7i) - log(fc 2 r? 2 ) = 




It is interesting to observe that the rate of spread in this case is three times slower 
than the rate of spread when the Richardson-Kolmogorov cascade dominates and one 
tries to collapse with outer variables. Hence, the spectral spread observed in figure [T4"k is 
in agreement with the view that regular grid turbulence at the Reynolds numbers of this 
figure obeys Richardson-Kolmogorov interscale dynamics. However, the high-frequency 
behaviours in figures [Mb . [T5b and [Tok fall within the uncertainty defined by (|A 8[) and 
we are therefore unable to conclude whether our fractal grid-generated turbulence obeys 
complete or incomplete self-similarity even though it is clear that L u « A is a good 
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approximation. By complete self-similarity we refer to the property that i<ii(fc,a:) = 
u' 2 L u f*(kL u ) is exact at all frequencies and by incomplete self-similarity we refer to the 
property that deviations to Fn(k,x) = u' 2 L u f*(kL u ) can appear at the very highest 
frequencies. We stress that this docs not imply Kolmogorov scaling even if these high 
frequencies may be collapsed by Kolmogorov inner variables for the simple reason that 
L u oc A and therefore e is not proportional to u' 3 /L u . 
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